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AN- EXPERIMENT TO EVALUATE SKYLAB EARTH RESOURCES 
SENSORS FOR DETECTION OF THE GULF .STREAM 

George A. Maul 
Howard R. Gordon 
Stephen R. Baig 
Michael McCaslin 
Roger J. DeVivo 

An experiment to evaluate the SIOOIAB Earth Resoutpes Package for 
observing ocean currents was performed in the Straits of. Florida in January 
1974. Data from the S-190 photographic facility, S-191 s^ctroradiometer, 
and the S-192 multispectral scanner, were coirpared' with surface observations 
made simultaneously by the R/V VIRGINIA KEY and the ilASA C-130 aircraft. . 

The anticyclonic edge of the Gulf Stream could be identified in the SKYLAB. 
S-190 A and B photographs, but the cyclonic edge was obscured by clouds.. The 
aircraft photographs were judged not useful for spectral analysis because • 
vignetting caused the blue/ green ratios to be dependent on the position 
in the photograph. The spectral measurement technique could not identify 
the anticyclonic front, but a mss of Florida ‘Bay water which was in the 
process of flowing into the Straits could be identified and classified. ■ No 
calibration was available for the S-191 infrared detector, so the goal of 
conparing the measurements with theoretical calculations was not accomplished. 
Monte Carlo simulations of the visible spectrum showed that the aerosol 
concentration could be estimated and a- correction technique was devised. 

'The "S-T92 scanner was not useful for detecting the anticyclonic front because 
the radiance resolution was inadequate. An objective cloud discrimination 
technique was developed; the resiYLts were applied to the several usefifL 
oceanographic channels to specify the radiance ranges required for an ocean 
tuned .visible multispectral scanner. 


1. INTRODUCTION 

An inportant problem in physical oceanography is the location of the 
boundaries of surface c’lrrents. Many techniques have been proposed to study 
such boundaries from space, but the actual process of extracting the correct 
infomation .from satellite data is in the early stages of development. 
SKYlAB, with its several types of sensors (NASA, 1975), afforded the oppo 2 ?- 
tunity to test three techniques simultaneously: photography, spectroscopy, 
and multispectral imagery. 

1.1 Background and Purpose 

Ifejor ocean currents are known to have several observable surface 
features that neke them disxinguishable from the surrounding waters. The 
Gulf Stream system is used as an exanple*to typify these changes because 
it is one of the most inportant ocean currents," and becaiise understanding 
its features permits one to study other current systems. 



Because of its subtropical origin, the current is typically wanner than 
surrounding waters and thus has a surface thermal signature that often can be 
detected ih infra 2 red imagery. The waters of the current are also much 
lower in biological productivity and hence there are fewer particles and 
biological pigment molecules in iiie Stream; this translates to water of a 
deep blue color. Conversely, the jxnctaposed water masses are frequently 
hi^er in biological productivity, and this can cause that water to be 
greener. Another feature of the current that naJces it visibly distinguish- 
able is caused by the large horizontal velocity shear, frequently the 
faster moving water in the current has a different sea' state than surrounding 
waters. Just as common are the many slick lines associated with the shear. 
Finally, the atmosphere above the Gulf S'tream is modified under certain 
conditions X'/hich can also give an indication of the current's location. 

There are several other features of the Gulf Stream that are potentially 
detectable by satellite altimetry and other microwave techniques . The 
measurement of these parameters will not be discussed in this report, the 
approach here being confined to visible and infrared wavelengths. The 
final goal of 'this experiment is to con'tribute to the multi-parameter deter- 
mination of the location of the Gulf Stream, of which visible and infrared 
measurements of radiance are an integral part. 

1.2 Test Site 

The experiment was conducted in the Straits of Florida along a sub- 
orbital- track in. -the vicinity -of Key -Vfest, Florida. In this channel, the 
Gulf Stream runs approximately perpendicular to the satellite ground track. 
This would maximize the change in oceanic variables while minimizing the 
impact on -the data acquisition facility onboard SKYLAB. Furthermore, the 
logistics of performing the sea truth from a 65 foot vessel in- January 
wea^ther made the choice of a semi-protected body of water, mandatory. 

l^drographic conditions in the test site are controlled by the 
location and intensity of the Gulf S'tream (also called, the Florida 
Current in 'this vicini-ty). The cyclonic edge, defined as 'the left hand 
edge facing downstream, has horizontal excursions of approximately 50 km; 
that is, at some times of -the year the current's edge my be found 20 km. 
south of Key West and at other times 70 km to the sou'th (kkul, 1975). The 
location of the cyclonic front determines the location' of the major 
hydrographic fea-hires of the S-traits, Materials from Florida Bay are -also 
known to flow into the straits and at times become entrained in the current. 
The detection of occurences involving mixing of Gulf Stream and Florida Bay 
waters is of fundamental inportance to the dispersal of na-tural and manr 
introduced mterials. 


1.3 Prior Investigations 

A general review of remote sensing of ocean color was given by Hanson 
(1972), and Maul (1975) discussed the application of visible spectroscopy 
to locating ocean current boundaries. ’ Gordon, in a series of papers (e.g. 
Gordon, 1373; Gordon and Brown, 1973; MaxH and Gordon, 1975) discussed the 
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spectra of upwelling irradiance as a function of the optical properties of 
the water as calculated by Monte Carlo simulations; those studies are 
directly related to the current boundary location problem because the 
spectrum of light changes from Gulf Stream water to coastal water. • Tech- 
niques for determining ocean chlorophyll (e.g. Baig and Yentschj 1969; 
Mueller, 1973; Duntley et al . , 1974) are also related to current boundary 
determination because pigment forming molecules, along with suspended 
materials, affecj: the light specxrum. 

Remote sensing of ocean currents in the infrared region of the electro- 
mgnetic spectrum has been attempted for many years (e.g.: Wamecke, et 
al. , 1971; Hanson, 1972; Richardson, et al . , 1973). Hosr?ever, there have 
been questions concerning the radiative transfer model dependency of the 
atmospheric correction (J^ul and Sidran, 1972; Anding and Kauth, 1972) that 
have awaited SKYLAB to be addressed. Adequate atmospheric correction 
techniques are required for ocean current boundary determination using once 
or twice daily observations because conpositing of images is required in 
order to fill in the areas covered by clouds; conposites must be based on a 
ccmmon measurement, that of the sea surface temperature itself - 

SKYLAB provided the first opportunity to evaluate photographic, 
spectrometric, and raultispectral imagery in a specific experiment designed 
for current boundary location. It will be seen that each instrument has ■ 
unique advantages, disadvantages, and limitations. It is the intent of 
this report to objectively evaluate each technique and provide recom- 
mendations for future equipment and measurements. 

2. GROUND TRUTH DATA 

This section gives the details of how the ocean surface data were 
obtained, calibrated, and analyzed. In many cases the surface optical 
measurements are useful indicators of the properties of the water that need 
to be. measured, . iMs is because the theory is well' ahead of llie m^sure- 
ments and adequate instruments are not yet designed or built. In the case 
of spiectrcHiieter measurements, the ideal observations are in fact physically 
impossible. 

2.1 Cruise Report 

The at-sea obseirvations were - designed to provide simultaneous measure- 
ments of the ocean while the aircraft and satellite transited the area. 

Since the speeds of the three vehicles are mismatched, the assumption 
must be made that the oceanic conditions are a steady- state for 12 hours or 
SO. While it is recognized that this is not strictly true, it is necessary 
in view of resources available. 

Weather conditions for the experiment were not ideal, with psartly 
cloudy s]<ys and moderate seas. Wind was from 045° at 4ms~^ and remained 
steady most of the day. Air temperature ranged fran 23.8°-26.1°C,wet bulb 
values were 23.0°C most of the day. Visibility was 20 km except in a rain 
shower at 1700 GMT where it dropped to less than 5 km. Barometer was steady 
at 1025 mb until 1900 GMT when it abruptly dropped to 1023 mb and 2 remained 
so thereafter. Wave hei^t was' 1/4 meter until 2100 GMT when it abruptly 
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increased to 1/2 meter; period remained 4 seconds throughout the day. 

Underway operations on the Virginia Key included gathering data on 
ocean salinity , chlorophyll-a concentration, surface nutrients, seawater 
scattering properties, sea sixrface temperature by bucket and by a 
continuous radiometric profile, and ocean temperature down to 450 m with 
, expendable bathythermographs; ship was lt>ve-to for spectrometry observa- 
tions. Collection of this data started at 24° 39'.1 N, 81° 08’.1 W at 
1253 GMT^ 8 January 1974. This point is about 7 kir . SSE of Marathon • 
in the Florida Keys. The track was directed SW and ended at 23° 33 ' .2 
N, 81° 55’ .5 W at 0150 GMT 9 January, This is 41 km . off the north 
coast of Cuba on the evening of the same day. 

The Virginia Key travelled at a speed of 8 knots while on track. 

The boat stopped only for spectronetry stations; all the trackline 
profile data was collected while making speed. 

• 2.2 Ttackline Profile Data 

Figure. 2.1 is a plot of the trackline data, after reduction. The 
stipled profile below the 22°C isotherm shows the bottom profile along 
the track. The arrow shows where the 22°C isotherm crosses a depth of 
100 m, a point interpreted to be approximately 15 km south of the 
boundary of the Gulf Stream (I^ul, 1975). The 2(45) curves show the ^ 
volume scattering function at 45 degrees, and are in units of meter 
steradlan (m"^ sr”^). 

A detailed discussion of the data collection, reduction, and 
interpretation- follows. 

a) Chlorophyll-a CCl-a) 

Cll-a concen-trations were obtained continuously by measuring 
the fluorescence of Cl-a when exposed to blue li^t. The data are 
• reported as if all the pigment forming molecules (including pheophytins) . 
were chlorophylls. The continuous record was obtained using a Turner 
fluoroneter, Model 111, which measured the fluorescence of surface water 
drawn' "thru a continuous-flcw intake system. This me’thod is as described 
in S-trickland and Parsons (1968), wi-th -the addition of a bubble -trap. 

In order to calibrate "the continuous record j three’ discrete samples were 
obtained by fil-tration and measured against a known standard after the 
cruise. This also was as described' by S-b?ickland and Paieons (1968) , 
and uses "the SCOR/UNESCO' equation. Table A.l in Appendix A shows the 
times- and positions of the "three samples . 

The large gap in the Cl-a -curve on Fig. 2.1 is due to a coiribina- 
tion of drift on a particularly long spectrometry station and a delay in 
turning on -the fluorcraeter after leaving "the station. The three other 
' breaks in the record represent a change in fluorescence- while stopped; ‘‘ 
for short spectrometry stations. The degree of variability in Cl-a 
concen"tration over the short distances indicated • in 'the. record shows "the 
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DEPTH (M) T (“C) SALINITY (%o) Cl-a{mgm*3) 


/S (45“) at 546 nm 


)S (45®) at 436 nm 


Bucket Temperature 


h/''v\ Radiometric Temperature 


22® C isotherm 


24® 39.1 N 
81“08.l'W 


23®33.2'N 
81®55.5'W 


Figure 2^1 Surface truth profiles across the Straits of Florida on 
8 January 1974. The profiles j from top to bottom are: Continuous 
chtorophyll-a (mg m'^; discrete salinity (°/qo); discrete volume 
scattering function at 45 (m~^sr~^); discrete thermometria temperature 

(°C)s continuous , radiometric temperature in 10.5 - Z.5.5pm band (°C)s 
discrete depth of 22°C isotherm. 
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desirability of a continuous record- instead of discrete sanples as a 
source of the profile. The hi^ values at the northern (left hand) end 
of the line occur over the reefs of the Florida Keys. 

In addition to the three discrete surface samples j discrete samples 
at various depths were obtained while stopped at two of the spectrometer 
stations . This „was achieved by acquiring water at the various depths 
for filtration and measurement later with the surface sanples. Times, 
depths -and positions are given in table A.l. 

b) ' Salinity (S°/oo ) 

A total of ten salinity samples were obtained while on the track- 
line. These were surface water samples which were bottled for 
measurement after the cruise. The times and positions of the salinity 
sanities are given in Appendix A, Table A. 2. The salinity profile in 
Fig. 2.1, which starts at 1200 GMT, is a strai^t-line plot of the ten 
•values obtained. 

c) Volume scattering C g ) 

The volume scattering function is a measure of the amount of 
scattering at various angles by a sample of seawater irradiated by a 
beam of light- Jn this case a single angle of M5° was measured, for a 
beam with a blue filter (436 nm) and a beam with a green filter (546 • 
.rim).., using -a .BriceriPhoenix li^t -scattering -photometer. g(45) was 
calculated, using : 

a(45°) = a TD D(45°) ' t SIN 45° 

TT h D(0°) 

there a is the ratio of the working standard diffuser to the reference 
standa3od diffuser, TO, is the transmittance of the reference standard, 
diffuser, h is the dimension of the irradiated element, D is the deflec- 
tion of the gaivonometer, and t is the -transmittance of -the neu-fcral 
density filters. 

Measurements were obtained by collecting water sanples wi-th PVC’ 
sampler bo-ttles. A' total of thirty surface samples -were measured, and-, 
at five stations samples were collected at various dep-ths. Values, 
times and positions appear in Appendix A, Table -A. 3. The (curve in Fig. 
2.1, which starts at 1200 GMT, is a s-trai^t-line plot of the surface 
values. 

d) Bucket Teirpera-tures (T^) 

As is customary, bucket tenperatures were acqLiired at each XBT 
cast. Additional bu<cket temperatures were acquired at spectroime-try 
stations and at samplings for scca-ttering measurements- The curve. in 
Fig. 2.1 , which starts at 1215 GMT, is a strai^t-line plot of the 28 
total tenperatures 'Obtained. See Appendix A Table A. 4 for values, times 
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and positions. 

e) Kadianetric Temperature 

A continuous sea surface temperature profile was obtained using a 
Bames, Model PRT-5, precision radiometric thermometer. This radiometer 
has a special 10.5-12.5 um filter that approximates those in the SKfLAB 
multispectral scanner. The -instrument’s voltage output was converted 
to temperature based on a calibration perfomed in March, 1974. .The 
profile in Fig. 2.1 begins at 1215 GMT. The large gap in the temperature 
profile matches the gap in the Cl-a profile and exists for the same 
re^on. 


■f.) Expendable Bathythermograph (XBT) 

A total of 23 XBT casts were made using the Sippican XBT system with 
450 m probes. These casts provided the information to plot the depth 
of the 22° C isotherm. The point where the 22°C isotherm crosses the 
100 m depth (arrow in Fig. 2.1) is taken to mark the zone of maximum 
horizontal velocity shear of the Gulf Stream. This happened at 
23° 40.9’ N, 81° 51.0’ W, about 57 3<m north of the coast of Cuba. See 
i^pendix A, Table A. 4 for times and positions. 

2 *^ <*M *JIJ 'I > piAtTJA 

■The spectral signature of the - ocean >and the sky at various -points 
along the' trackline was measured in the visible range of li^t. The 
vessel was stopped for the spectrometry stations; six stations were 
occupied for a total of 25 spectra. Table 2.1 shows the times and- 
positions of these stations. A detailed discussion of the instrument, 
data: reduction- and wavelength calibration follows. 

a) Ins-trument 

The ins-trument was a Gamna Scientific, Model 2400 SR, spectro- - 
radiometer in a special watertight case.- The Model 2400 SR scans in a 
wavelength range of 350 to 750 nanometers (nm) by rotating a high efficiency 
diffraction grating which faces a narrow aperture slit. It has a wave- 
length, accuracy of _+2.5 nm. For -this experiment the^ instrument was set 
to scan from 375 to 725 nm with a Wra-tt^ 2b filter installed over the 
entry slit. This filter effectively cuts transmission, below 400' nm, 
insuring -that the results would not contain secondary diffraction return. 

An opal glass diffuser plate was used as a cosine (Lambertian) collector; 
all -measurements are irradiances. The data were recorded on a dual 
cihannel strip chart recorder. 

b) Data Reduction 

Fig. 2.2 is an example of a -typical spectral scan. It is a 
scan of ocean upwelling li^t from the station at 2210 GMT (see Table 
2,1). The dashed curve with the many peaks is -the original unsmoothed- 
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TABLE 2.1 


TimeCGMT) 


Position 


Spectrometty Stations 


Solar Zen.ith 
'Angle 


c 

nj 


0) 

o 

o 


n3 

L 

O -P 
• O 
• OJ 
O Cu 
■ -I2: to 


1406 

24° 

38.9 


81° 

08,0 

W 

64 ® 

3 

1 

1603 

24° 

30.7 

N 

81° 

16.3 

W 

51° 

3 

2 

1815 

24° 

19 .”1 

N 

81° 

27,3 

W 

47° 

3 

2 

1956 

24° 

08.1 

N 

81° 

34,6 

w 

56° 

4 

1 

2ioo 

24° 

04,0 

N 

81° 

37.0 

w 

65 °- 

2 

0 

2210 

23° 

58.8 

N 

81° 

40.4 

w 

76° 

3 
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WAVELENGTH (nml 


Fd-gure 2.2 Example of upwelZdng speotval dvradiance h&fove (broken tdne) 
and after (soldd tdne) fdttering^ to etdmdnate the effect of ocean 
surface gtdtter vardatdons due to surface waves. ' : 




data. The peaks, are due to the ci^aiiges in the angle of the water surface 
relative to the instrument as waves pass underneath. This imposes a 
fairly regular periodic variation over the general trend of the spectral 
return; all of the spectra acquired on this cruise showed this varia- 
tion to some extent. It is necessary to remove these peaks if the data 
are to he useful. 

Digital low-pass time series filtering techniques were used to ^ve 
the spectra the smooth appearance shown by the solid curve in Fig. 2.2. 
These techniques, to be e3q)lained in some detail below, permit an ob- 
jective filtering of the imwanted periodicities with a minimum, loss of 
si^uLfioant data trends. The filtering was performed on a UNI VAC 1108' 
using tte FESTSA (Herman and Jacobson, 1975 ') software system at the 
Atlantic Oceanogmphic and Meteorological Laboratories. 

The spectra were digitized off the strip chart at an interval of 20 
points per inch, which was sufficiently fine to retain the shape of the 
original traces. For the trace in Fig. 2.2, this provided 248 data 
points. As different wavelength drive speeds were used on different 
stations, this number varied from scan’ to scan. 

In order to choose a filter that successfully res^ves high frequency 
energy with a minimuin loss of significant trends, it is .important to 
identify the periods of the high frequency signals. Fig. 2.3 contains a 
plot (light 5 broken line) of the relative strengths of. various periodicities 
.of the, .spepfipmetry .scan shown ,in .Fig. 2.2. This is a pow^ spectrum of 
the spectrometry scan using TUKEY's method (see Herman and Jacobson, 

1975).' One can see strong periodicities at approximately 6,7,9,20 and 
30 data points per cycle. These are the dominant hi^ frequency signals 
which give the original scan in Fig. 2.2 its sawtooth appearance. 

-The response of the filter chosen to remove- these hi^ frequency 
signals is shown by the heavy, solid curve in Fig. 2-3. . This response 
is in terms of a ratio of the contribution various frequencies make to 
the form of the original trace, to the contribution the same frequencies 
are allowed to make to the form of the filtered result. Thus in the 
example chosen, no contribution is allowed for periods smaller than 
about 9 data points per cycle (10 nm), fiill contribution is allowed for 
periods greater than about 90 points per cycle (100 nm), and half contri- . 
butioiis are allowed at about 25 points per. cycle (30 nm), which’ was the 
longest period of the major peak in Fig. 2.3. 

]ji general,, all filters were chosen, to remove the short period 
(hi^ frequency) signals in the same way. Power spectra of different 
spectrometry scans did not always closely resemble, eachi other however, and 
each filter had to be chosen on the basis of an individual inspection of 
eacdi scan. 

c). Wavelength Calibration ^ 

Irradiance and wavelength are- indicated by separate voltage _ 
outputs. The wavelength voltage is produced by a potentiometer directly 


- 10 - • 



PC^WER (H*. ARBITRARY UNITS) 



SAMPLE INTERVAL 


^ spectrum of upwelUng speatral ivvadicmaa shorn in 

motion omssd hy 

st^faae waves reneot^ng spsoularly; the heavy line is the response 
of the Fovrver fvltev to low-pass the data. 
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connected to ’the diffraction grating » voltage varying with angle. . In 
addition, an irscribed wavelengtfi scale is connected directly to the 
potenticmeter. Thus it is possible to canpare the voltage. of the wave- 
length output as recorded by whatever strip cj-iart recorder is used with 
the wavelength indicated by the scale. Such a conparison was perfornied 
on the strip chart recorder used during the cruise, and is the basis of 
the wavelength calibration employed for this data. The drift free 
nature of the grating-scale design permits this type of calibration. 

Recorder outputs were ccmpared with scale readings at 5 nanometer 
intervals over the entire' wavelength range used in this experiment. 

A third-order polynomial was fitted to the resulting numbers to obtain 
an equation giving wavelength in terms of the position within the spectra 
on the time axis of the strip chart:. This form of equation was chosen 
because it does not require a fixed number of data points for all 
spectra; the only ipequirement is that the digitizing interval ronain 
constant. 

CciTiparison of this calibration with 3<nown spectral lines; indicates 
the calibration is within 3 nanometers of true wavelength over the whole 
range of visible li^t. 

3. PHOTOGRAPHY 

The use of color photographs to measure phytoplankton concentra- 
^J.ipns .in ..riatural, waters ,is based on, the .argument that the photographic 
material responds in a quantitative, reprxxiUcible manner to the variance 
in the li^t field of the vjater. The variance is associated in a fixed 
manner with the concentration of the phytoplankton,- its distribution 
with depth, and with its species and nutritive history. For- a number of 
years there have been applications of these variance techniques to 
remote sensing of the ocean (e.g. Baig and Yeutsch,-. 1969; Mueller, 1973), 
The SKYLAB experiments provided the opportunity to extend some of the 
techniques developed from laboratory tanks and low-level aircraft 
flints to synoptic mesoscale coverage. The field. program (section 2.1)* 
was to provide surface truth for the variance analysis ; the- SKYLAB 
photography was to provide the photograpMc products . 

3.1 Measurements 

.A photograph is merely the record of the integral of the intensity 
of the illuminant falling on a sxjbject and the reflectivity of the 
subject. In a color-photograph a third variable is introduced in the 
spectral properties of both the illimninant and the subject. A color 
photograph is satisfactory if it produces, ,in a viewer.' s eye, a response., 
similar to that produced by the actual' subject. The satisfactory spectral 
response brought about by the color photograph is a' result of the eye ’ s 
inability to distinguish between a pure spectral .source of li^t'and a 
mixture of such sources. ' A color photograph does not produce in each 
pixel (picture element) the exact spectral reflectivity of the corres- 
ponding spot of the subject (with the exception of a subject which ’is 
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BLUE/6REEN RATIO 



Figicee 3, t Exampte of the effect' of vignett-ing in aircraft film data 
on the percent trdnsirfittance of blue (450 nm) and green (550 Tim) 
Zight^^ and on the hZue/green ratio, AericCL ooZor positive fiZm was 
used in the ECS camera. 
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itse!Lf a color photograph) . Instead the photograph produces in each pixel 
a mixture of three colors which the .eye perceives as a single color, 
and 'it produces only these three colors: this is called a "raetameric" 
match. 


A color can be thou^t of as a vector in n-space, with pure light 
as ^ the origin. «An infinity of coordinate systems nay be created around 
this vector, but once one of the coordinate axes is fixed the others 
also become fixed. Coiimon varieties of color films need only three 
colors to reproduce the variety of colors seen in the real world. To 
the eye these colors look like yellow, magenta, and cyan (blue-green). 
These colors are the coordinate axes of the color space. Every other 
color is an unique combination of these three primary colors. If the 
subject is conposed of two different colors then the eye will perceive 
it as if it were a single color. The eye in this case performs 
the inetameric natch. A color photograph will 'do 'the 'same: If the color 

of a subject is changing then the change my be noted as differing 
quantities of the three dyes in a color photograph of the subject. 


The utility of monitoring the changes in dye concentration of a 
color photograph can be carried a step further. It has been shown that, 
provided the continuous spectrum of the subject has previously been 
measured, the dye concentrations can be used to generate a new spectrum — 

T*T“l i-i” C'n-vn-i t v ti rmi 


T nco N 

/ * 


spectrum must have been part of a ’’training set” of spectra for which 
the color photographs exist, or must be any combination of the original 
training set spectra. Then, throu^ a multivariate analysis and regres- • 
sion technique, a synthetic spectrum can be generated using only the dye 
concentrations. The technique is especially useful when the concen- 
tration of one of .the conponents of a mixture is changing. Tank, (Baig 
and Atwell 1975) and low-level aircraft flints (Baig, 1973) have 
amply demonstrated that phytoplankton concentration in natural' waters 
can be easily and accurately measured with the 'technique. If the spectra 
of the phytoplanktors are. not of immediate interest Ihen "the multivariate 
reduction is not necessary. The problem then reduces to a correlation 
between the concentration of phytoplankton in a training sample and "the 
variation in the dye concentrations in “the color photograph. 


3.2 Data Analysis 


The first photographs to be analysed were 1he 9" color 'transparencd.es 
frcra "the aircraft serial cameras. On a nuniber of frames 'there were sub- 
jective color differences. However, similar differences were -noted on 
frames in which •the color of -the subject area would have been expected 
to be rjniform. A densitcme'tric analysis of one of "these frames revealed 
a variable blue/green ratio as a function of radial distance from the 
principal point. These data are presented graphically in Fig. 3.1. 
A'ttenpts to use "these data to- ’’correct" data from other frames of aircraft 
transparencies were not successful. This vigne'tting problem was so 
severe -that the images were displayed on -the non-linear portion of ihe 
film’s D-log e curve. This in'troduced an unknown non-linear error which 
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could not be corrected without an in-fli^t calibration with targets of 
known spectral response/ 

Because the scale of the satellite photos is so much siraller than 
that of the aircraft photos, discontinuities such as fronts and eddies 
are recorded in only a small area of the photo. By conparison, similar 
elements have t® be recorded over large areas of single aircraft photos, 
or even in a sequence of such aircraft photos. The practical effect 
of this scale difference is that variations in film density related 'to 
position on the photo can be ignored where the area of interest only 
covers a. small area of the photo. Of course comparisons between areas 
widely scattered over the photo are still subject to the problem of 
spatial density in the photo. 


All of the pertinent SL-4 duplicate films were analysed on a hybrid 
transraissometer. The li^t table and associated apertures, filters, and 
diffuse acceptor are from a Welch Denischron densitometer. The li^t 
sensor and associated electronics are a Gamma Digital Photometer. The 
transmission of each of the three color filters and the white li^t 
setting was calibrated with a non- silver standard step wedge traceable 
to N.B.S. standards. Sucdi a step-wedge is a better approximation to the 
actual attenuation characteristics of color film than the usual silver 


grain step wedges. This is because color films do not have any silver 
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of the illurainant. Data in the following paragraph are reported as the 
percent fraction of ‘transmitted ii^t rather than as density, as trans- 
mission ratios will have a more meaningfijl interpretation 'than density 
ratios. Standard deviations (a) follow each ratio. 


The first area analysed was a plume of water off the coast of Cuba, 
in frame 97, roll 64 SL-4, near reseau (Fig. 3.2). The average 
Blue/Green transmission (B/G) ratio is 5.9, a +_ 0.1. Just to the left of 
the interface of this plume with the Gulf Stream water the B/G ratio 
changes to 7.0,cr + 0.1: in the Gulf Stream water- In frame 98 a similar 
plume on the FlorTda Keys side of the Gulf Stream has a B/G ratio of 
6.0, cTjf 0.1, while the Gulf Stream VJater immediately to the right of the 
plume has a- B/G ratio of 6.8, cr+_ 0.1 (Fig. 3.3). This particular plume 
shews in three frames, 97, 98, and 99. The B/G ratios for the plume are 
5.5, 6,0, and 6.2, respectively, the B/G ratios for the Gulf Stream 
water adjacent to the plume are 6.3, 6.8 and 7.1, respectively. Thus, 
while the absolute values of the ratio are changing, the difference 
between the plume ratio and the Gulf Stream water -ratio is nearly constant 
fxon frame to frame. 

Both filtered panchromatic films SO-022 showed some apparent density 
changes in the same areas as those in rell 64. Roll 65, filtered to 
pass 0.6 to 0.7 micron li^t showed a transmission change from 40.0 x 1(?^ 
to 34.0 X 10”^ on going frem the plume off Cuba to adjacent Gulf Stream 
water. Roll 66, filtered for. 0.5 to 0.6 micron li^t showed no trans- 
mission change between these two areas , both noted as 34 x 10 , 
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EEPEODUCmiLITY OF THE 
PEIGINAL PAGE IS POOR 



Fi-aure 3.2 S-190B -panGhromatio (SO-022) photograph of the Straits of 

Florida and the northern coast of Cuba. Fhe box in the lomr left 
brackets the antiayolonic front of the Gulf Stream j and deft-nes the 
area where densiometrio measurements were made. 
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Figure 3.3 S-V90B panahromatic (SO-022) photographj of the Straits 
'' of Florida and the western Florida Keys. The box in the upper left 
brackets a plume of water from Florida. Bay^ and defines the area 
where densiometrio measurements were made. 
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Neither of the two b/w I.R. fihns showed any -transmission differences 
among the areas analyzed. The color I.R. film did show some differences 
which are considered ro be statistically significant. The plume off of 
Cuba showed a B/G ratio of 7.0, while adjacent Gulf Stream water showed 
a ratio of 7. -5. ^ The plume off the Florida Keys showed a B/G ratio of 
8.5, while acros's the interface of the Gulf Stream water showed a ratio 
of 7.2.- It should be noted that the B/G ratio of the color infrared 
film is really a ratio of the green to visible red. radiation. 

The 'conclusions which can be drawn fron 'this limited data set are 
■that a significant variation in ocean color* can be observed by changes - 
in dye concentrations in color photographs of the scene. By considering 
"the surface 'truth the evidence tends to favor the variation in suspended 
cihlorophyll as the most probable cause of -the color variation. 

3.3 Disoiission 

Tt is immediately apparent from the data ■that ■the transmission 
ratio technique is a useful means of analyzing variations in color of 
satellite-derived photography. At the same time the data might have 
been more useful had certain precautions been taken. Reference .is 
specifically made to 'the aircraft-derived photography. ' To achieve a 
flat spectral response across the film, the associated optics should 
have been fitted with anti- vignetting filters. The space craft cameras ' 
-suffered to -a lesser extent wi-th the same problem-. -In -this case each of . 
■the associated optics had been calibrated so that "the error in 'transmission 
was known. There is however, no indication "that such care was taken in 
preparation of the subsequent duplicate images. V/hile care was taken to 
ensure that duplicate grey scales were reproduced. 'at the same levels as 
those on the on-board films , apparently no account, was taken of -the 
variation in illumination across .the print head of 'the printer. ‘All of 
■these problems taken together substantially reduce 'the possible intei>- 
conparisons -that might have been a-ttempted. 

The photographs in Figs. 3.2 and 3.3 are bo-th S-190B products which 
have been enhanced by printing on hi^ contrast film. Exposure levels 
were set to saturate "the details in "the non-oceahic features . This is a 
•trial and error technique that extracts narkedly more low radiance level 
information. The change in texture marking differing sea states in Fig. 

3.2 is not measureable using the densitometer technique, but it is 
clearly noticeable to "the eye. The boundary be'tween' 'the two levels of ' 
radiance is probably -the anticyclonic edge of -the Florida Current. The ' 
detection of the features in Figs. 3.2 and 3.3 by "the S-192 scanner is 
discussed in section 5. - • 

■4, SRECrROMETER EKPERIKENT’' 

The SKTLAB S- 191 steerable spectroradiometer was to be used in- 'this 
experiment to study changes in the visible C 0.4-0. 7 y m) and infrared 
(7.0- 14.0 ym) spectra' of the ocean across -the current’s cyclonic boundary. 
The plan, was for the crew to acquire a cloud-free oceanic -area with the 
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S-191 looking 45° forward of nadir, and track that site mtil 0°. 
Thereafter, the spectroradiometer was to be locked into nadir viewing 
across- the cyclonic fronr and up into the waters- of Florida Bay. 

The experiment proved to be not very successful for several reasons: 
while the crew did as the plan said,. the data acquisition camera (DAC) 
was turned off and the exact tracking data (angles, times, locations) ■ 
were never recorded; the calibration of Hie S-191 infrared detectors is 
not known; the visible region data radiances values do not agree with 
theoretical or observed values reported by other investigators. The DAC 
problem will, be discussed in this introduction, and -the infrared and 
visible data in sections 4.1 and 4.2, respectively, below. 

Location of the data is made difficult by -the lack of DAC ou-tput. 

The voice log was the best clue to what actually was done by the crew. 
From the transcript of the voice tape, at 16:29:33 GMT the pilot had the 
S-191 set at 35° looking forward along the track, although all crew 
instructions were to set the S-191 target acquisition at 45°. The word 
"thirty five" was nob clearly audible however, and the pilot may have 
followed the instructions sent up to SKYLAB just prior to 'the pass. At 
16:29:35 Gf^, the pilot reported tracking a clear area of water; the 
assigned start time was 16:29:33. The exact time of reaching nadir is 
difficult to tell from the voice log, however 16:30:45. is the approximate 
tiji^. 


The location and time of the nadir point was calculated from geo- 
*metricai considerations assuming a spherical earth with a radius of 6378 
km and a satellite altitude of 443 km. If the nadir angle was 45° at 
16:29:33 GMT, the position of the point tracked was 23° 53' -4 N, 81°58*.0W 
the time of arrival of the spacecraft over this point from the best 
available positioning data is 16:30:41.1 GMT, which is in good agreement 
with the voice log estimate. The message sent up to the crew had the 
finish of the tracking at 16:31:05. 

The S-192 line s-trai^tened data show that the position given above 
is in the middle of a clear ocean area and it appears reasonable to have' 
tracked this as. the site. Ihe vehicle was over Florida Bay at 16:30:55 
according.-to the. S-132, but the pilot commented at 16:31:10 that they 
were going across -the jfeys. This discrepancy camot be accounted for 
unless the Florida Ifeys were observed well after 1be spacecraft transit. 

If the above analysis is correct, then according to -fche ground 
truth data in Fig. 2.1, the Gulf Stream was probably never acquired by 
the S-191. The position of the 22°C isotherm at 100 meters dep-fdi" ■ 
indicator was 25 km SW of the point where the pilot tracked a clear 
area. Although the exact location of the front cannot be identified in 
the ship track data it appears that it was also SW of the nadir tracking 
point. Maul (1975) reported the mean separation between the indicator 
isQ-therm and the front to be 11 km in this area and that data further 
supports the contention that the S-191 did not obtain spectra in the 
Gulf S-tream. The objective of analyzing the change in spectra 
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across- the front cannot be accomplished with these data. 

4.1 Infrared Eadiance 

The infrared experiment was designed to study the accuracy of 
atmospheric transmission models. This objective cannot be accomplished 
because the calibration of the infraired detector is an unknown function 
of wavelength (Barnett, personal communication; Anding, 1375). Several 
relative tests were made however, which provide sohb informtion on the- 
atnpspheric transmission model dependency, and these will be discussed 
below. 


4.1.1 Theoretical Calculations 


Ehiitted infrared radiation (7 yiTi <1 ^14 ym) leaving the earth 
passes throu^ the atmosphere before detection at the S-131 sensor. The 
atmosphere modifies the infrared radiation by absorption and, to a very 
minor degree, by scattering. Details of the theory aipe given by Chandrasekar 
(1960) and recent reviews on its application to oceanography were given 
by Hansen (1972) and Maul (1973). The radiative transfer equation 
throu^ an absorbing but non-scattering atmosphere is: 


N(9,x) = e(e',X) L (T ,X) t(9,X) 


L(Ta,p ,x) 3T(p,e,A) p 


f 

Jo 

+ p (6', X) N (9",X) t(6,X) 


(4.1) 


•3p 


as 

<■ * 

where 6,9’ ,X" are. the nadir angle, angle of reflectance, and angle of 
incidence, respectively. Radiance (N) at the satellite is wavelength 
dependent, and is a function, of the surface blackbody radiance (L) , the 
emissivity of the surface (^) and the transmittance of the atmosphere 
(t); these three parameters describe the absorption of emitted blackbody 
mdiation by the atmosphere. The second term in the equation, the 
integral term, describes the atmospheric (a) modification of the radiance 
as a- function of pressure (P) . The third terra describes the contribution 
of the reflected ( p) atmospheric radiance at the surface , again as 
modified by transmittance. 

The theoretical calculations discussed herein are an extention of 
the model used by Maul and Sidran (1973) which uses the transmissivity 
data of Davis and Viezee (1964). The area of interest-, is the 10.5 - 
12.5 wn band which is used on many spacecraft, including the SKYIAB S-192 
multispectral imager. In this spectral interval e> 0.99 at low nadir • 
angles; hence pC=l-e ) is very small and equation 4.1 irey be written 


N(0) = ({»(X) L (Ts,X) t(9,X)dX 




, 4i(X) L (Ta, p,X) 3£(p,9, X) dpdX 

3p 


(4.2) 
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ATMOSPHERIC PRESSURE (mb) 


RELATIVE HUMIDITY (%) 



AIR TEMPERATURE (“0 


Figure 4,t Vertical profiles of atmospheric pressure and relative 
humidity taken at the time-^ of SFILAB transit. The dotted lines 
on the relative humidity profile are the cloud-free estimate of 
atmospheric moisture. 
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The filter function ((j> ) is zero outside the interval discussed above. 

The radiance may be converted to equivalent blad<body temperature by 
inverting the Planck equation CD which has been integrated over the 
same 10.5 < *<12. 5 interval. 

The calculations were carried out on the AOML computer. A special 
radiosonde was released by the Key West office of the National Weather 
Service at the time of SKYIAB transit (see Fig. 4.1). Before computing 
the radiative transfer from" a radiosonde, the data must be inspected to 
insure that no clouds are in thr. path of ascent, in order to compute a 
cloud- free radiance. Clouds are readily identified by their charac- 
teristically high- relative humidity and isothermal temperature. There is 
evidence of clouds in the data in Fig. 4.1, and calculations were 
na.de" to test the effect of clouds . 

IWo cloud layers are in evidence, one centered at 744 nib and one 
centered at 671 mb. Clouds are characterized by a sudden increase in 
relative humidty and a snail (near zero) lapse rate. An equivalent 
clear sky estimate is made by assuming the clouds are absent; ^ the 
estimated relative humidity profile in the clouds region is- ^ven by 
-the dotted curve. The calculated equivalent blackbody tempera"tures for 

T 298.15°1C ara: 
s 

Wavelength Observed (Appendix B) , Cloud- free Equivalent 

ilm 293. ‘22° K 293. 85° X 

12.5m 290.46° X . 291. 28° x 


The differences in this case are snail, 0.63°x at 11pm, and 0.82°x 
intergrated over the 10. 5-12. 5pm region where "the S-192, -NOAA-4/5 and 
SMS-1/ 2 observe. Other experience with this type of cloud-free equiva- 
lence has been as high as 5°x over the Gulf of Mexico. 

4. -1.2 Comparison of S-191 and Models 

As stated in section 1. 3, the wavelengths chosen for. the two - 
channel technique, (Anding and Kauth, 1970) of atmospheric correction 
depend on the radiative -transfer model. SKYLAB was to be used to study 
•that question but since -the calibration of -the S-191 infrared detector 
is unknown, -the problem cannot be investigated. 

The mean sea surface temperature along "the -trackline was 25.0°C. 

This value has been used in the calculations shown in Fig". 4.2. The 
Davis and Viezee (1964) model does not include absorption due to--the 
ozone miolecules which show up as a maximum at 9 . Sum" in -the S-191- ob- 
servation. The conparison shows -that ozone- does not affect -the 10.5- 
. 12.5Pm window and hence is not a factor in the S-192 infra 2 ?ed scanner 
data. At ll.Oym, -the apparent difference be-tween Jhe observed and cal- ■ 
culated equivalent blackbody temperature (Tgg) is--3.5°C. . This seems to 
be -the approximate error estimate of o-ther SKYLAB:' investigators (personal 
c»nmunicatiohs), but no conclusions can be. drawn. 
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Figure 4.2 Spectral infrared radiance observed by the S--191 spectra-^ 
radiometer (dashed line) ^ and calculated by the ozone excluding 
model of Davis and Viezee (fine solid line), Eeavy solid lines- 
are blaakbody curves. 
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Since "the calibration uncertainty is wavelength dependent, the data 
at 11. till were studied to determine the shape of the nadir angle depen- 
dence curve. In the lower half of Fig. 4.3 is a least-squares fourth 
order polynomial fit to all the observed radiances as a function of time 
(dashed line). „ Since the same ocean spot was to be tracked, radiance 
should be a function of nadir angle up to 16:30:45 GMT. The mximum on 
this curve (arrow) is at 16:30:19. The upper curve is the theoretical 
calculation using equation 4.2 for the same atmosphere in Fig. 4.1 and 
for T =25®C. Nadir angles were confuted using a start time of 16:29:35 
GMT (45°) and a stop time of 16:30:45 GMT (0°) following the discussion 
in section 4.0. The match in the curves maxima would be approximately 
coincident if the tracking started with a 37° nadir angle, which is in 
agreement with the voice tape transcript. 

The S-191 spectroradiometer uses a series of detectors which cover 
a segment of the spectrum. In some regions these overlap and ambiguity 
often exists as to which detector to use. In the sections to follow, 
those detectors (channels) chosen were as recommended in the NASA reports • 
on instrument performance. It is suggested that only that data which is 
well calibrated be reported to non-instrument engineering investigators 
in the future. 

Bartlett (personal communication) cautioned against the use of 
radiometer' channel A1 in the S-191 (see again Fig. 4.3) Accordingly a 
second fourth-order polynomial (solid line) was fi-tted to the channel A6 
data only. The rms spread of the radiance about this polynomial is 
4.48 vW cm^ sr^ This correspxDnds to a noise equivalent temperature 

difference (NEAT) of 0.3°K at 289.65°K. Since the a-tmospheric attenuatiqn 
tends to dirninish surface gradients , this results in a calculated NE A T 
of approxiirsately +_ 0.6°C in. Tg for this model at 11 ym on this day. 

The equivalent blackbody temperature at the maximum in the polynomial is 
16.5°C at the top of the atmosphere. This implies a tempera'ture correction 
of 8.5°C which is not unreasonable ‘for a tropical winter atmosphere 
whose precipitable water vapor is 3.6 cm (cf Maul and Sidran, 1973). 

4.2. VISIBLE RADIANCE - 

The visible experiment was designed to study the accuracy with 
vdiich spectral changes across oceanic fronts can be observed and inter- 
preted from satellite altitudes. Unfortunately, the strongest front 
expected in the exp>eriment area was missed by the S-191 so this objective} 
as discussed before, cannot be accomplished. Hcwever, during the course 
of this, work, a theoretical technique for recovering the "ocean -color 
spectrum" through the atmosphere was developed. This is discussed in 
detail below and an attempt is made to conpare the predictions of the 
idieory with the S-191 data and the associated groimd truth. 

•4.2.1 Theoretical Calculations 

It is clear that "the full potential of oceanic 2 ?emote sensing from • 
space in the visible portions of the spectrum can be realized only if 
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Figiwe 4,3 Radian< 2 & at ttym as a funat-ion of nadiv angle of the same 
ooean spot as tracked on the S-l9l (data), Channel A1 data to the 
left of the vertical dashed line; channel AS to the right. Dotted ■ 
line is fourth order polynomial fit to all data; fine solid line is 
fourth order polynomial fit to AS data' only. Heavy solid line is 
the calculated nadir angle dependence. 
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■the radiance which reaches the top of 'the a’tjnosphere can be related to 
the optical properties of the ocean. To effect this, the radiative 
•transfer equation must be solved for the ocean-a'tmosphere system with 
collinated flux incident at the top of the atmosphere. In such calcula- 
tions 'the optical properties of the ocean which must be varied are the 
sca'ttering phase function (P© ( 9 )) and •the single scattering albedo 
Cwq; defined as” the ratio of the sca’ttering coefficient to ’the total 
a’ttenuation coefficient). Furthermore, unless the ocean is assumed to 
be homogeneous, the influence of vertical s’truc'fcure in •thes.e properties 
must be considered. To describe “the cloud free atmosphere, the optical 
properties of "the aerosols and their variation with wavelength and 
alti’tude as well as the ozone concentration' must be known. Considering 
•the ocean for ’the present to be homogenous, 'the radiance at the satellite 
can be related to ’the ocean’s properties by choosing an a-tmospheric 
model and solving the transfer equation for several oceanic phase functions 
and cuo^s at each wavelength of interest. The number of separate computa- 
tional cases required is then ’the product of the number of phase functions, 

' the number of values of Wq, and the number of wavelengths . Even if ’the 
multi-phase Monte Carlo me’thod (MPMC, Gordon and Brown, 1975) is used, 
the Wq resolution of Gordon and Brown (1973) would require a number of 
simulations equal to ten times the -number of wavelengths for each 
atmospheric model considered. It is possible, however, to obtain the 
necessary information wi-thout modeling the ocean's optical properties 
in such detail. 

The model is based on an observation evident in results of 
coiputations given by Plass and Ka.’ttawar (1969) and by Kattawar' and 
Plass (1972) on radiative -transfer in the ocean atmosphere system, 
namely, that when -the solar zeni-th angle is small, ' -the upselling radiance 
just beneath the- sea surface is approximately uniform, (i.e. ,-not s-trongly 
depenaent on viewing angle) and hence determined by the upwelling irradiance. 

This observation is utilized 'in simulations of oceanic^remote sensing 
situations by assuming that a fraction R of -the downwelling photons are 
absorbed. The ocean is then -treated as if there is a lambertian reflecting 
surface of albedo R just beaneath the sea surface. In this case Gordon 
and Brown (1374) have shown •that any radiome-tric quanti'ty Q can be • 
wri-h+en 

■ Q R 

. Q = Q + 2 (4.3) 

1 1-rR 

Qj is "the con-tribution to Q from photons which never pene-trate 
•the sea surface (but may be specularly reflected from -the surface). Q 
is .the contribution to Q from photons which interact with the hypo- 
•thetical "labertian surface" once for the case R-l. r is "the ratio of 
the number of photons interacting with -the "lambertian surface" twice, 
to. -the number of photons interacting once, again for R=l. 

Using equation 4.3 any radiome-tric quanti’ty can "then be ccaiputed as 
a function of -R. Physically -the quantity R is -the ratio of upwelling to downwelling 
irradiance jiist benea'th -the sea surface and is known as -the reflectance 
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Table 4.1 , The Three Ocean Scattering Phase Functions 


e 

(deg) 

KA 

CxlO^) 

KB 

CxlO^) 

KC 

CxlO^) 

0 

10924 

10171 

9521 

1 

,4916 

4577 

4285 

5 . 

573.5 ■ 

534.-0 

499.9 

10 

169.3 

15 7.7 

147.6 

20 

29.5. 

29.39 

29 . 31 

30 

12.56 

11.95 

11.42 

45 

3.059 

3.661 

4.189 

60 

1.092 

1.577 

1.999 

75 

0.546 . 

0.915 

1.190 

30 

0 . 344 

0.661 

0.952 

105 

0.311 

0.641 

0.928 

120 

0.317 

0. 7.3 2 

1.094 

'135 

'0.410 

0.829 

1.309 

.150 

0.492 

1.017 

1.618 

165 

0.579 

1.261 

1,856 

180 

0.617 

1.357 

. 1.939 
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function' [R(0,-)'3 in the ocean optics literature CPreisendorfer, 1961). 
Spectral measurements of the reflectance function R{X ) have been presented 
for various oceanic areas by IVler and Smith (1970). Henceforth, R(x ) 
will be referred to as the "ocean color spectrum." 

■A series of Monte Carlo computations have been carried out to see 
if an approximate simulation (ASl), using this assumption of uniform 
upwelling radiance beneath the sea sxorface, yields results which agree 
with cojiputations carried out using an exact simulation (ES) , in which 
the photons are accurately followed in the ocean as well as the atmos- 
phere. The Monte Carlo codes used in Gordon and- Brown (1973, 1974) were 
modified by the addition of an atmosphere. -The atmosphere consisted of 50 
layers and includes the effects of aerosols, osone,- and Rayleigh scattering, 
using data taken from the work of Elterman (1968). The aerosol 
scattering phase functions were computed by Eraser (NASA-GSFC, personal , 
communication) from JiLe theory assuming an index of refraction of 1 . 5 
and Deirmendjian's (1964) "haze C" size distribution. Also, to determine - 
the extent to which the vertical structure of the atmosphere influences 
the approximate simulation, a second approximate simulation (AS2) was 
carried out in wliich the atmosphere was considered to be homogeneous; 
i.e., the aerosol scattering, Rayleigh scattering, and ozone absorption 
were independent of altitude. The oceanic phase functions in the ES are 
based on Kullenberg's (1968) observations in the Sargasso Sea, and are 
gxven Table 4.1. 

■KA is roughly an average of Kullenberg ’ s phase function at 632.8 nm 
and 655 nm, and KC is his' phase function at 460 nm. KB is an average of ■ 
KA and KC. These phase functions show considerably less scattering at 
very snail angles ( 6 <I°)than observed by Petzold (1972) in other clear 
water areas; however, the exact form of the oceanic phase function is 
not very inportant, since it has been shown (Gordon, 1973) to influence 
the diffuse reflectance and R(0,-)only through the back-scattering 
probability ( B) ^ 

B = / p (9 ) sin 0d0. 

In all of the conputations reported here, the solar beam incident on the 
top of the atmospshere, is from the zenith and with unit flux. At visible 
wavelengths, the variable atmospheric constituent which will most- strongly 
influence the radiance at the top of the atmosphere is the aerosol 
concentration, so the conputations have all been carried out as a function 
of the aerosol concentration. 

Table 4.2 gives a sample of the results comparing the upward flux 
at the top of the atmosphere at 400nm in the three simulation models 
(ES, ASl, and AS2) as a function of the aerosol concentration. N, 3 x N, 
and 10 x N respectively refer to aerosol concentrations in each 
layer of 1, 3, and 10 times the normal concentration given by Elteman. 

400 nm is <^osen because in the visible portion of the spectrum it is 
the wavelength at which the atmospheric effects are expected to be most 
severe. The ES case uses w =0.8 and phase function KC. The values 



Table 4,2; Comparison of the flux at the 
top of the atmosphere for the 
ESa ASls and AS 2 simulations. 


Aerosol 

Concentration 

ES 


, ASl 


AS 2 


N 

0-222 

(+.002) • 

0.224 

(+.001) 

0.226 

(+.001) 

3xN 

0.274 

(+.003) 

0.273 

(+.001) 

0.275 

(+.001) 

lOxN 

0.423 

(+.004) 

0.426 

(+.002) 

0*425 

( + .002) 
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of R used to effect the AC computations were taken from 1±ie EC compu- 
tation of this quantity, however if R is taken .from • 

R = 0.0001 + 0.3244X + .1425x^ + .1308x^ (4.4) 

idiere x = u^B/ (L~Wo(l-B) which, according to Gordon, Brown and Jacobs. 
(1975), reproduces the in--water reflection function for. the corres- 
ponding case but with no atmosphere present, the results of the AS 
computations agree with those listed to within 0.2%. The numbers in the 
par^thesis next to each flux value represents the statistical error in 
the flux based on the actual number of photons collected in each case. 

It is seen that ES and AS simulations generally agree to within the 
accuracy of the computations. Notice also the excellent agreement 
between the Asl and As 2 fluxes. 

In Tig. 4A the comparison between 'the. ES, ASl and AS2 upward radiances 
at the top of the atmosphere is presented for the three aerosol models. 

The step-like curve in the figure is for ES, the solid circles for ASl, 
and the open circles for AS2, and Vpis the cosine of the angle between 
Idle 'nadir and the direction toward which 'the sensor is viewing. The 
radiances in Fig. 4.4 for the ES cases are accurate to about 3% in the 
range y= 1 to about 0.4, while for the AS cases the accuracy is about 

1 Si T.TT+Vinfj .. — i^— 

«-4.Vrfis.> y L4 u:v.^v^ I <4 lO 

again agree for all the aerosol concentrations except within the range 
iP. .about J3„..3.; i._e.., viewing near the -horizon. It -is. felt that these 
computations demonstrate that 1die transfer of 'the' ’ocean color spectrum 
through the a'tmosphere can be studied with either the' ASl or AS2 nicdel 
as long as radiances close to the horizon are not of interest- Further- 
more, from the reciprocity principle (Chandrasekhar, 1960) -^e nadir 
radiance, when the solar beam makes an angle So- with the zeni-idf,' can 
be found by multiplying the radiance I( 4 ) in Fig. 4.4 by v vhere y is 
taken to be cos ®q. This implies that as long as the sun is not. too 
near the horizon, 'the ASl and AS2 methods of, computation can be used to 
determine the nadic radiance at the tip of the a'tmosphere as a function 
of the ocean’s properties through eq_uation 4.4. The fact that the AT2 
model (honpgeneous atmosphere) yields accurate radiances is very important 
in remote sensing since it inplies that only 'the total concentration (or , 
equivalently the total optical thickness) of the aerosol need be determined 
to recover the ocean color spectrum from satellite spectral radiometric 
data. ' ' 


It should be pointed out that "these results also strongly suggest 
■that R( X ) is the quanti-ty relating to the subsurface conditions , which 
can be quantitatively detemnined from space, and hence, is "the most 
na"tural definition of the "ocean color spectrum". Moreover, it has been 
shown (Gordon, Brown, and Jacobs, 1975) that R( X ) is not a strong 
function of the solar zenith angle (the maximum variation in R(0',-) with 0^^ 
is of the order of 15% for 0 ■^9o-i60°)in con'trast wi-th o’ther definitions 
(Curran, 1972; Mueller, 1973). 
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Figure 4,4 CompartsoTL hetween ES (siep Zihe curve) ASl (soUd ctvpZes) ond 
ASS (open airales) upward: radiances at the top of the atmosphere for 
an ocean with Wq ^0,8 and phase function KC and an atmosphere with . • 
a normat (IxE)^ three times normal (3xN) and ten..,times normal (ZQxN) ’ 
aerosol concentration. . * 



The only way spacecraft data can be used to obtain inforration 
concerning subsurface conditions (such as concentrations of chlorophyll, - 
suspended sedinents, etc.) is through determination of R(X ). [It is 
assumed here that the relationship between RC X ) and the ocean constituents 
is well ]cnown, while in fact much work still renains to be carried out 
before such a relationship can be established ]. This can be effected 
by applying equation 4.3 to radiance I(v ) at the top of the atmosphere 
with the sun at the zenith, which yields, 

R ICii) 

Um) = I Cu) + ^ 

1 1 - rR 


II (y ) and I ( y ) are presented in Figs. 4.5 and 4.6 respectively for 
the three aerosol models discussed above as well as an aerosol free 
model (QxN) and a model with seven times the normal aerosol concentration 
(7xN). For the eases considered, PsO.5 and since R is usually about 0.03 
to 0,10 at this wavelength vje can rewrite equation 4.5 approximately 
as 


i(y)= I, (y) + Ri^(y-) ■ . (4.5’) 

so R 0 ^ Ky) ~ T (y) C4.6) 

“ I Oy) i . • 

2 

implying the reciprocity theorem to equation 4.6 it is found for nadir 
viewing, that 

-R=^ ^nadir - (4.7) 

where 4o is the cosine of the solar zenith angle Qq ’ Noting again that R 
is = 0.10 it is seen that the difference between and yo-i(4o) 

must be small which implies that the accuracy in R will be limited by 
knowledge of li(yo). Since 1^(4^) depends strongly on the aerosol con- 
centration, it is absolutely necessary to be able to determine the 
aerosol concentration if an accurate value of R is desired. Curran (1972) 
has suggested that this can be accomplished by observing the ocean 
(assumed free of white caps) in the near infrared where R(^ )= 0. In 
section 4.22 of this report, Curran’s suggestion is utilized to find 
the aerosol concentration from the S-191 data. 

Before trying to apply the relationships developed here to the S-191 
spectra, there are several important irrplications of the theory which 
shall be discussed. Noting that I^( 4 ) and I 2 ( y ) depend only on the 
direction of the incident solar beam, the properties of the atmosphere 
and ocean surface, but not R, if it is assumed these latter properties 
remain essentially constant over horizont^ distances large compared to 
those over which R changes significantly, one can directly relate changes 
in I(y ) to changes in R. From equation 4.5’ 
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Figure 4.S I^(v-) as a function of y for various aerosol concentrations 
■wavelength of calculations is 400 nanometers. 
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3lCy) I (y) • 

— m 2 

Figure 446 shows that 81/ 8 R is not an extremely strong fmction of the aero- 
sol &oncentration for concentrations up to three times normal , and viewing 
angles up to 35? from nadir. This suggests that horizontal gradients in R 
can be estimated without knowing the aerosol optical thickness with great 
accuracy. 

Using equation 4.5 to relate changes in radianceCAl(y) ) to changes 
in R( AR)‘, i.e. 

I(y) = C3I( y)/9R]ARc:I (y)AR. (4.8) 

— 2 

Equation 4 . 8 a enables a determination of minimum radiance change the 
sensor must be able to detect for a given ^R. For exanple, suppose 
■that observing at y= 0.85 it is desired to detect a 5% change in R for 
clear ocean water at 400 nm CR~.l) through an atmosphere with 3 times 
the nomal aerosol concentration. Figure 4.6 shows that !I^(0.85) is 
about O.ll and noting that the extra terrestrial flux 400nm is about 
140 W cm"'^ nm* , we find from equation 4.6 that AK0.85) is 0.077 Wcm'^nm^sr 
In a similar way radiance changes can be related to A R for a nadir-vievring, 
sensor and any solar zenith angle. As mentioned previously frem xhe 
reciprocity principle, 

I _ = I(yo) 
nadir 

where y^= cos^o? ®o is the solar zenith angle, and KPo) is the radiance 
at the top of the atmosphere seen by a sensor viewing at ® q when the sun 
is at the zerdtii. Following through with the same arguments that led to 
equation 4.8a it is found that 

Al =y CalCy J/3R) Cyo>AR. ' (4.8b) 

nadir 0 0 ° j 

Clearly, for a given AR, decreases ‘substantially with increasing 

solar zenith angle due to the presence of the y^ factor in equation 4.8b. 

For example, with a 3 times normal aerosol concentration, a nadir 
viewing, sensor would have to have about 2.5 times more sensitivity at 
Qq = 60° as conpared to 0 to detect the same R. 

The above examples indicate how the theory (ASl) can be used in 
the design of a satellite sensor system for estineting some ocean property 
such as the concentration of suspended sediments or organic material. 
Specifically, one must first determine the effect of the property to be 
investigated on R, then based on the sensitivity desired, find AR, 
and finally, use equation 4.8a or b to find the rainimum radiance change 
•the sensor must be capable of detecting. If the sensor has a limited 
dynamic range, then equation 4. .5 can be used with equations 4.8 a or b 
to aid in the sensor performance design trade-offs. [Unfortunately 
at this time, relationships be-tween R(A) and sea water consti'tuents 
are not well established. 3 ■ . . 



i I i t 1 i 1 


i ■ I i 



"HAZE C" 

HENYEY- 

GREENSTEIN 


..,^© 0 = 0.6 

"g^07~~ 

^08 


I I I I ( I I J ( t 1 1 I I I I t 


0 20 40 60 80 100 120 140 160 180 


0 (degrees) 


Figure 4. 7 Comparison between the "haze C" and various Henyey-Greenstein 
phase functions characterized by asymmetry parameters 0. 6^ 0, 7 y and 
Q,8y as a function of scattering angle (Q), ' 



Considering the fact 1±iat we have only used the "haze C" aeixDsol 
phase function (which is clearly only approximately dharacteristic of 
the actual aerosol scattering) it is natural to inquire how strongly 
the coirputations of I (y) and 1 2^1^^ presented in Fig. 4.5 and 4.6 depend 
bn the shape of the aerosol phase function. To effect a qualitative 
understanding of the influence of the aerosol phase .function, conputa- 
tions of I and have been carried out using the well known Henyey- 
Greenstein^(HG) phase function 

» 

. (0) = a-g^)/4TT 

Cl+g^-2g cos0)'^/^ 

where the asymmetiy parameter g is defined ' according to 

'' g = 2 / P(0) cos 9 sin 0 d0, 

J o 

and 0 is the scattering angle. Since g for the '*haze C" phase function 
is 0.690, conputations have been made withPf^Q(0) for g values of 0.6, 

0.7, and 0.8. Figure; 4.7 ccmpares these Ppjgfe )’s with the "haze C" 
phase function. The H6 phase function for g = 0.7 clearly fits the 
"haze C" phase function quite well in the range 5° <0 <140°; however, as 
is well known, the HG formula is incapable of reproducing phase 
functions computed from Mie theory in the extreme forward and backward 
directions. The KG phase functions wixh asynmetry parameter 0.6 and 
0.8 are seen to he substantially different from the I'haze C" distri- 
bution at nearly all scattering angles. On the basis of Fig. 4.7 it 
should be expected that I and I ^ coirputed with P^qC® ) will be in 
close agreement with the "haze C" computations only for g close to 0.7. 
Figures 4.8 and 4.9, which compare the results of computations of I,- and I 
respectively for %g(9) with g= 0.6, 0.7, 0.8 (step like lines) and 
idle "haze C" phase function (solid circles) for the normal aerosol 
concentration, show that this is indeed the case. It is seen that except 
for apparent statistical fluctuations , the HG phase- function for g = 0 . 7 
yields values of 1 1 and I 2 in good agreement with the "haze C" coirputa- 
tions. This suggests that the detailed structure of the phase f met ion 
is not of primary importance in determining and 1 2 , and it may be 
sufficient for remote sensing purposes to parameterize the phase 
function ty g. 

To get a feeling for the importance of variations in the phase 
function in the remote sensing of ocean color, consider the effect of 
changing the aerosol phase function frem an HG with 'g=0.6tog=0.8 
over an ocean with R = 0.1. From Fig. 4.8 and 4.9 it is found that the 
normalized radiance at y = 0.85 (the assumed observation angle) decreases 
by 4.9 X 10 ^ , this decrease in radiance would be interpreted under the 
assunption of no atmospheric change as a decrease in R from 0.10 to 
0.056. This clearly indicates then that variations in the aerosol phase 
function in the horizontal direction could be erroneously interpreted as 
horizontal variations in the optical properties of, the ocem. It is 
however probably unlikely that except in extreme cases, the clear 
atmosphere iceanic aerosol phase function will exhibit variations as large 
as considered in t±ds example. Assuming that the aerosol concentration 
of the atmosphere can be determined, the uncertainty in the aerosol 
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Figure 4,8 Comparison hetneen computed for the ^'haze C" and Henyey- 

GreeHstein phase functions for an atmosphere liith a normal aerosol 
concentration, as a function of cosine B( v); ‘wavelength of calculations 
is 400 nanometers. ■ . • - 








(001 X 




phase function will still of course provide a limit to the accuracy 
with which the ocean colo 2 ?-spectrum can be retrieved from satellite 
radiance measurements. 

In summary then, the theory (ASl) leads to the natural definition of 

[R(0,-) as a function of wavelength] 7 as the "ocean color spectrum". 

The determination of subsurface oceanic properties from space can thus 
be divided into two problems; 1) the determination of R(A) from satellite 
radiance measurements, and 2) the establishment of relationships betv?een 
R(x) and the desired ocean properties. Since the method of coTiputation 
conveniently separates the radiance into a coirponent that interacts with 
the ocean (I 2 ) and a component due to reflection from the atmosphere and 
sea ^surface CI^, it is easy to relate changes in radiance to changes' R(X) . 
It is found that for vievTing angles up to 35° from nadir, I 2 is a rel- ' 
atively weak function of the aerosol concentration for concentrations up 
to 3 times normal. This suggests that spatial gradients of ROt ) can 
be determined with only a rou^ estimate of the aerosol concentration. It 
is further found that variations in the aerosol phase function can strongly 
influence the interpretation of the radiance at the satellite. Clearly 
then, it is vital to understand -the magnitude of aerosol phase function 
variations . 

4-. 2. 2 Technique for Atmospheric Ctorrection 

As discussed in section 4.2.1 it is necessary to know the aerosol 
concentration in order to recover the "ocean color spectrum" R(^ ) from 
the nadir radiance .spectrum observed at the .satellite. In this section 
a method based on Curran's sijggestion of using the near infrared radiance 
to determine the concentration is developed and applied to the S-191 
data. 


The method involves finding a band of wavelengths in the near 
infrared for which 'the absorption by ozone and water vapor is negligible. 
Since the Rayleigh scattering by air is very small -in the near infrared 
the greatest contributor to the optical thickness at the wavelength in 
question is the aerosol. It is found that at 780 nm ozone and water 
vapor do not absorb significantly, and the Raylei^ scattering only 
contributes about 0.023 to the total optical thickness of the atmosphere, 
implying that aerosols play the dominant role in the radiative transfer 
here with the normal aerosol concentration yielding an optical thickness 
of about 0.2. Also since R(0,-) for wavelengths greater^than about. 

700 nm is essentially zero, the upward radiance at the top of the atmosphere 
at 780 nm simply becomes 


1(4) - 


= 


vhere 


(mW cm"^ 4'ni^ 


F^is the solar irradiance 
mosphere and is taken from Kondrat'ev. 


), at the top of the at- 
(1973). Using the the reciprocity 
principle, for nadir viewing and any solar zenith angle (9 0 ^ 9 ^ 
can be written nadir 


I 

nadir 




I j(vi) and have been computed for aerosol concentrations 0xN,lxN 



2xN and 3xN at: 400 nm, 500 ran, 600 ran and 780 ran and the results for the 
0'^ and liiN computations are presented in Appendix C. Using I^Cu)' ; 
for 780 ran and noting that the S-191 nadir radiance was recorded with 0^= 
40° the upward radiance at the top of the atmosphere for nadir viexving is 
found to be 0.323 and- 0.938 mWoiT^ sr"^ for aerosol concentrations of 

QxN and IxN res;^ectively. Using the S-191 radiances at 780 nm for' 
the nadir viewing spectra taken on Jan. 8, 1975 at 16:30:45,75 GMT 
(spectrum A) and 16:30:52.2 Gid? C spectrum B) which respectively were' 

0.64 and 0.72 ir£V cm^ sr“^wn^ , it is found that theory suggests the 
aerosol concentration was 0.51 xN for spectrum A and 0.64 xN for spectrum B. 
In order to compute R(A ) from the S-191 data the assun^tion is made that 
the variation of the aerosol extinction coefficient with x^7avelength is 
exactly as given by Elterman. 

4.2.3 Recovery of R(X) from the S-191 data 

As mentioned above, in order to recover R(X ) from the S-191 data 
it is necessary to assume that the variation of the aerosol extinction" 
coefficient with wavelength is identical to. that .given by Elterman. Also, 
since Ij (4) and I (u) at 780 nm were derived using' the *’haze C” phase 
:^ction for the aerosols, the assumption this phase function is correct 
is implicitly being made. With these assumptions the nadir radiance 
at the top of the atmosphere has been computed at 400, 500, 600 and 730 
■ran for -aerosol concentration OxN -and IxN assuming that R(x) is. zero, and 
these 'are presented in Table 4.3 along with, the radiances from spectra A 
and B. 


Since the aerosol concentration is known to be between OxN and IxM 
it appeal^ that the' S-191 data at 400 ran is in error. It is virtually . 
impossible for the nadir radiance to be less than that for a OxN atmos- 
phere. (It should be noted that the discrepancy here is not small, 
i.e. the S-191 radiances at 400 nm appear to be too small by more than 
a factor 'of 2.) The radiances at the other wavelengths listed in 
Table 4,3 seem to be reasonable and these will be used to -estitnate 
R'(^) using equation 4.7. The results of this are shown in Table 4.4. 

’ ■ \ 

It is seen that R(X ) is negative except in the spectral region 500- 
550 ran where the values shown compare well with the Tyler and Smith Gulf 
Stream data for R(0,-). As discussed above, the '40 0_ ran data is appar- 
ently in error,' however the data at other wavelengths appear realistic,. 
so the negative R(0,-) values are probably due to 1iie assumptions 
■that the "hazeC" phase function characterizes the aerosol, and 'that 
the spectral variation of the aerosol scattering coefficient is correctly 
described by Elterman’ s data. It should be clear' "that considerably more, 
experimental work is required in order to test the ability of the theo^ 
discussed in 4.2.1 to accurately find R(^) from the satellite radiance. 



Table 4,3 


Wavelength 



I 

nadir 

Tr^oI 


(nm.) 

ImW cnf-^pm 


mw cm '^: 

pm sr ^ 



Q 

OxN 

IxN 

Spect A 

Spect B 

400 

157 

5.61 

■ 6.58 

2.30 

' 2,50 

500 

201 , 

2.93 

4*10 

4.13 

4.39 

600 

184 

1.25 

2.13 

1^55 

1.67 

780 

125 

0.323 

0,938 

0.64 

0.72 



Table 4.4 


W5.VS XsTi^*fch 

R(0 , 


Crun; 

Spectrum A 

Spectrum B 

400 

-0.274 

, -0.268 

450 

-0.0242 

-0.0235 

500 

0.0318 

0.0370 

550 

0.0295 

0.0303 

600 

-0.00715 

-0.00547 

780 

0 

0 
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5. MULTISPECTRAL SCANNER EXPERIMENT 


SKYLAB’s imiltispectral scanner was a unique design that had 13 spectral 
channels of data' spread over the visible and infrared bands. The system 
used a conical scan which had the advantage of keeping the atmospheric 
path length the same at all times. The visible region of the spectrum 
-(0.4 - 0.75yra) was divided into 6 channels, each about 0.05 nin wide. 

IVo reflected infrared (0.75 - l.Oym) channels and one in the emitted 
infrared (10.2 - 12.5 ym) were also provided. Hie channels useful to 
oceanography, which were examined in this experiment, are summarized in 
table 5.1. 

LANDSAT-1 has been shown to have several useful applications of 
•visible region imagery to marine science (Maul, T9 74). The much finer 
spectral resolution of the S-192 provided an opportunity to expand those 
results to ocean current boundary determination and to test if the 
lower wavelength ( 0.5 ym) intervals were useful throu^ the intervening 
atmosphere. 


5,1 S-192 Data 

~Srl92 idata were collected from 16:29:22 GMT (over the open sea just 
north of the Ctiban coastline) to 16:31:04 GMT (over the mainland Florida . 
coast north of Florida Bay). All channels listed in table 5.2 were 
carefully examined in the analog format provided by 'NASA to the principle 
investigator. 'The data in the images were compared' with the S-190A 
and S-190 B photographs to see if what is interpreted in section 3.2 as 
the anticyclonic edge of the current could be detected. This feature 
was not observable in the standard data product. 

The cyclonic edge of the stream appears to be obscured by clouds. 

This is often a useful means of locating the edge of the current but 
unfortunately made the objective of directly sensing the edge an im- 
possibility. 

An unexpected opportunity to evaluate the S-132 however, is evident 
in the photographic detection (section 3) of a mass of water from Florida 
Bay flowing south into the Straits of Florida just west of Key West. 

This water is observed to be milky in appearance and somewhat greener an 
color. No ocean surface spectra were observed inside or outside of the 
plume of Florida Bay water, although it could have been easily accomplished • 
if the SKYIAB crew had observed the feature and radioed the observation, so the 
the Ship knew of its presence. Upwelling spectral irradiance reported 
by Maul and Gordon (1975) probably describe the essential features of the 
plume and water in the straits. 

An intensive effort was made by Morris’ (NASA-JSC) , Johnson (Lockheed- 
JSC) , and Maul (NOAA-AOML) to extract the two features discussed above 
using the computer enhancement facilities at. NASA-JSC. After approximately 
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TABLE 5 .,1 


BAND 

DESCRIPTION 

RANGE 

(lim) 

1 

Violet 

0.41 - 

0.46 

2 

Violet-Blue 

0.46 - 

0.51 

3 

Blue-Sreen 

• 0.52 - 

0.56 

4 

Green-Yellow 

0.56 - 

■ 0 ..Bl 

5 

Orange-Red 

0.62 - 

0.67 

6 

Red 

0.68 - 

0.76 

7 

Reflected Infrared 

0.78 - 

0.88 

8 

Reflected Infrared 

0.98- 

1.03 

13 


10.2 - 

12.5 
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Figure 6. Z Fewer spectrum of the radiance in the unfittered S-Z92 

conical format. Significant noise is noted every . 5 , lO-$^ and 

6 data po^nts. • . 
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10 hours, of nechine -time on both conical and line straightened data, the 
feature described as the anticyclonic edge could not be identified, 
althou^ it is clearly brought out in the photographic enhancements 
(see Fig. 3. 2 ). Further effort to bring out the anticyclonic edge 
vas judged to be tmwarr anted and attention was turned to the plume 
feature which is visible in Fig. 3. 3, and which preliminary enhance- 
ment at the NASA-JSC facility showed to be a useful area in which to 
worko 


Prior to developing a general computer, enhancement technique, 
the data were examined for periodic features in a TUKEf spectrum. 

Figure 5.1 is a spectrum of data specially provided for this experimenr 
that was to be high pass filtered only; the calibration of the S-192 is 
considered a high pass filter. Significant periods at about 15 
data sample intervals is noted in these conical data as' has been 
reported (Schell, Philco-JSC, personal communication 1975.) The line 
strai^tened data (see Fig. 5.2) have been band pass filtered to remove 
this 15 data sample periodicity . The wavy patterns near the edge of 
clouds is idle result of filter ringing. 

5.2 Ctonputer Enhancement 

Computer enhancement of S-192 data was an objective of the experi- 
ment. The technique described below is a step toward autonEitic 
detection ox clouds in multispectral data. The goal is to use a near 
infrared channel (channel 8 in this case) to specify where cloud free 
areas are for analysis of sea surface temperature, or ocean color. 

Channel 8 ( 0.98-1. 03 y m) is selected as the cloud discrimination 
channel because there is a maximum in the atmospheric transmissivity at 
this wavelength, and amxijiium in the absorption coefficient of water. 

The hi^ absorption coefficient of water at ly m causes the -ocean- 
surf ace to have a very low radiance when compared with land or clouds. 

Thus there should be two modes in the frequency distribution of radiance : 
one mode for the clear ocean and another mode for land and/or cloutfe. 

An exanple' of such a bimodai distribution is given by the histogram in 
Fig. 5.4. 

In this figure, the low ocean radiances are clustered at the mode 
centered at N - 0.2 yW om'^ sr‘^ . The other mode, centered at 
N = 5.7 yWcm^sr”^ is a contribution of the clouds (there is no land in ■ 
this example) . If these modes can be identified and separated, a statistical 
identification of cloud- free ocean pixels can be mde. 

Cox and Hunk (1954) observed that the radiance reflected frcm -t±ie 
ocean is essentially Gaussian in character. The problem then is to fit 
a curve- of the form 

y = exp [-(N - N)^/2cr23 
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(5.1) 





Figure 5,2 $-Z92 tine strdightenedj filtered, scanner data over tUe Strarits 
of ’ Florida near the western Florida Keys., The S-192 channel -number is at 
the top of each appropriate panel. 



NORMALIZED FREQUENCY OF OCCURRENCE 



0 1 2 3 4 5 6 7 8 9 10 II 

RADIANCE {/xW cm-^sr"') 


Figicps 6. 3 Ristogvam (novmdti-zed to unity) of the vadiance in ahanneZ 8 
ovcT the seme ar-ea as figure 5.2 the primary peak at the left is clear 
oaeanj the broad peak centered at ? cnr'^sr'‘^is due to clouds 'and land. 
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to the data at the lower valued mode. In this equation, the nornal 
frequency curve Cy) is a function of iiie total number of observations (n) , 
the class interval (i), and the standard deviation Co); the pverbaron 
the dependent variable (N) denotes ensemble average. Fitting equation 
5.1 to the data is done in an iteration scheme that uses the lower valued 
mode as a first 'estimate of N (only the values 2a from the original _ 
ensemble are used in this first iteration; this eliminates much of the 
cloud contaminated data). After the first fit using a prede scribed N, 
the Sterne is to iterate the data using only ^2a of each new fit- When 
a (or N) changes less than 0.1% between iterations, the fit is considered 
acceptable and the cloud free pixels are defined as those between 
O^y*^ ‘*'3. This guarantees that 99% of the values around the lower mode 
are accepted and included in the ocean data. ■ Experience with this 
algorithm suggests that only 5 or 6 iterations are usually necessary 
for the scheme to converge. 

In Fig. 5.4, the data from the left hand portion of Fig. 5.3 
are plotted along with the fitted frequency distribution givoi by equation 
5,1. This fit required 5 iterations . Cloud-free ^ata are conservatively 
identified as all those whole N< 0.. 5 yW cm2 gpi (M+3a). Data from any other 
channel can now be statistically examined on a pixel by pixel comparison 
vn.th channel 8; only ”cloud-free” values go into the statistics. 

jCloud .free .data from .any other channel are analyzed for their mean 
and standard deviation. These calculations automtically provide the ’ 
limits over which 'the ocean data are to be stretched. Following the 
technique suggested by Haul, Chamell, and (Jualset (1374), .the formu- 
lation used by Maul (1975) is used here- The stretch variable ( C ) for 
a negative image is defined by: 

C = 0 for N >N-+Ka 

? = M C(N+ Kg) -M] for (N - icn )< N <(N+ko^ 

2ko 

C = M.for N <N- 

In this formulation M is the maximum value allowed by the Digital-to- - 
analog (D/A) output device, N is the mean radiance of the cloud free 
data, and is a constant. Considering N a continuous variable setting 
k= 2 would stretch 95% of the cloud free data over the full range of the 
photographic ehhancen^t device. ’ • 

The technique described above allows an objective specification of 
both the range of settings for ,an optimum enhancement of an ocean scene-, 
and a statement of the radiance range required ox an ocean color sensor 
under these conditions. In Fig. 5.6, the data from the sediment plume 
flowing Through -the Florida Keys are presented for channels 1-6,8, and 
13. These data are not line-strai^tened and they are only hi^ pass 
filtered. While the data are distorted geographically, they represent 
the best radiometric infometion from the S-192. 
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NORMALIZED FREQUENCY OF OCCURRENCE 



RADIANCE (/aW cttr® sr“' ) 


in Figure S.3. 

Todicmce dit^tulitn. ^ ‘Vprc^ation to the oUervfd . 
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figure. 5 . 5 Coniaai S-192 imagery of the same area as shown in figure 5. 2. 
CkanneZ 8 in this figure is a binary mask with alear ocean hlaak^ and' 
otovdSf tand, and other unwanted -pixels white. 















TABLE 5,2 


Channel 

Mean ± a 

Oceanic ' 
Hadiance Range 
mW cm^sn^Uirf^ 

1 

3,56 ± 0.45 

2.21-4.91 

■ 2 

4,71 ± 0.97 

1.80-7.62 

3 

4.07 ± 1.33 

■ 0.08-8,06 

4 

2.83 ± 1.46 

0.00-7.21 

5 

1.73 ± 0.92 

0.00-4.49 ■ 

6 

0,93 ± 0.43 

0.00-2.22 

7 

0.57 ± 0,29 

0.00-1.46 

8 

0.37 ± 0.22 

0.00-1.03 

13, 

0.80 ,± .0,03 

0.71-0.89 
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5 . 3 Discussion 


V7ithout actual spectra as surface truth it is not possible to adequately 
interpret the scanner data in the vicinity of the plume. Other spectral 
data (see Section 4.2 and ffeul, 1975) suggest that there should be a 
detectable difference in the data from channels 1 throu^ 5. Little 
information, let alone information difference, is contained 'in channel 1 
(see again Figs. 5.2 £ 5.6), These findings are in -agreeinent with Hovis 
(NASA-GSFC, personal corrmunication, 1975) who found that little or no 
infornELtion is detectable throu^ the atmosphere for wavelengths much less 
than 0.45 ym. The Kaylei^ scattering at shorter wavelengths is so intense 
that multispectral scanners such as the S-192 or the coastal zone color 
scanner destined for NIMBUS~G should not have a violet (0.41 - 0.46ym) 
channel. This means that the wavelengths which contain much of the oceano- 
■ graphic information are of little value at orbital altitudes. 


Analysis of the conical data used in Fig. 5.6 allows a specifi- 
cation of the radiance range encountered. Since an ocean color sensor 
should be allowed to satijrate over land or clouds, the range appears 
narrow. The range is based on + 3 cr , which describes 99% of the data 


Slice ui 


1^0. CUU I 1 I l^±. 


I. 






in order to expand the quantization commensurate with an acceptable data 
flow rate. ' Table 5.2 lists the radiance ranges. 


This range of values applies only to this data set, which represents 
low latitude, winter conditions. A similar- analysis of S-192 data from 
other 32 ?eas of the oceans and at other times would lead to an objective 
statement of the radiance range specification for an oceanic color 
sensor. Note that the lower value on the shorter wavelength channels is 
non-zero this ref leers the radiance of the atmosphere only. 


Neither the line-straightened (Fig. 5.2) nor tiie conical (Fig. 5.5) 
imagery capable of detecting any ocen thermal variations. The 

thermal data (channel 13) were processed in both negative and positive 
format in an attempt to show some of the 2°C range recorded in the 
sijrface observations (Fig. 2.1). In canparing the data, note ^ that the 
area covered is not identical due to geographical distortion in the 
line-strai^tened processing. 

The actual mask generated from using equations 5.1 is illustrated 
in the channel 8 data in Fig. 5.2 . Here the values 0 or 3 are assigned 
for contaminated or ocean (respectively) radiances- The channel 8 data 
in Fig. 5.5 were sinply stretched over the + 2 cr range by equations 5,2. 
A significant error in the automatic stretch units of other channels 
would occur if the extra step of 0 or 1 assignment -were not -mde , because 
of the exclusion of some ocean radiances. Note aiab that the effect of 
filter ringing is eliminated in the Fig. 5.2 nasking; this can also lead 
to wrong stretch limits if care is not exercised in application of the 
technique. Filtering should always be done after the data are masked. 
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The conical scan -technique is an unusual approach to multi spectml 
scarinihg. ' Data quality of the S-192 is judged to be poorer -than experience 
with the LANDSAT MSS which uses a linear scan. If the poorer S-192 data 
is inherent in the design of all conical scanners then it cannot be 
reconmended for the NIMBUS-G^ coastal zone color scanner- If this is not 
the case there is no reason based on this investigation co not consider 
such a future design. 


6. SUMMARY 

The objectives of -this experiment were to obtain simultaneous ship, 
aircraft, and spacecraft data across the Gulf Stream in the Straits of 
Florida in order to evaluate several techniques for remote sensing of 
this ocean current. Calibration of "the S-191 infrared detectors was not 
known; this precluded comparing -the a-tmospheric transmission models » 
which was an objective. Detection of "the current with the photographs 
was possible but the S-192 scanner was not useful in this respect be- 
cause the gain settings were not adequate for ocean radiances. A technique 
to determine atmospheric corrections for visible radiometers based on 
theoretical considerations was shown to be promising. Details of -(hose ■ - 
results axe enumerated below: 

1. Obsexvarion and data reducxion techniques for surface truth 
vessels include .measurement of ocean color spectra from 3 meters above 
■the surface. Objective filtering techniques were developed to remove 
periodic specular return caused by wave facets, (section 2.3) 

2. The limited photographic (S-190) data set in this experiment 
provides a baseline against which satellite photos of mere biologically 
productive waters may be cempaxed. Inter comparisons with aircraft 
derived photography may be more difficult. This problem was not re- 
solved in -this experiment because of unforseen variation in the exposure 
level of these photos. (section 3.) 

Non-uniformity in the satellite space-derived -photography is probably 
mostly due to -the effects of variable lens -transmission. It is recommended 
that the production of multiple generation photos ("dupes") be docu- 
mented in a manner sinilar to that of the ori'ginaT films. This should 
include some documentation on the printer light variation across -the 
platen, suspected to be a major source of the variability measured in 
the dupes used in this ejqjeriment. (section3.2) 

3. The data acquisition camera on SKYLAB was not -turned on during 
the experiment; this made any quantitative comparisons wi-th spec-tra 
iinpossible. There should be a positive interlock on all future missions 
to prevent a reoccurrence. Spectra were not obtained across "the Gulf 
Stream front because a) no data just prior to -the mission were obtained, 
and b) no direct ccniiiunication link from "the ship to the spacecraft 

could be set up. Fu-ture experiments must include a mechanism to communicate 
to surface -truth investigators -the position of variable fea-tures such as 


- 56 - 



cx^ean currents, (section 4.0) 

,4. Unfortunately the S-191 visible near-IR data could not be used 
to study the observability of oceanic fronts from satellite altitudes 
because the Gulf Stream front was missed. However, a theoretical method 
for recovering the "Ocean Color Spectrum" throu^ the atmosphere was 
developed. This was used to try and retrieve the "ocean color spectrum" 
from the nadir- viewing S-191 data, with limited success. The results 
agreed well with measuranents (Tyler and Smith, 1370) for wavelengths 
greater than 500 nm, but in the blue the S-191 radiance was substantially 
smaller than previous aircraft observations,' and even a factor of 2 less 
than theoretical predictions for an aerosol free atmosphere. This is 
either due to a mn .function of the sensor in -the blue or to 'the presence 
of a very s-trongly absorbing aerosol. Our present knowledge of the optical 
properties of marine aerosols does not appear to be sufficiently complete 
to effect a quantitative retrieval of "ocean color spectrum" from spacecraft 
data. Clearly further research on this problem is indicated. (Section 4.2) 

5. Ocean fea^tures that were visible in photographs (section 3.2) 
were not visible in 1he S-192 multispectral scanner data due to the lack 
of radiance resolution. ^The range of radiances needed to adequately 
observe ocean features is presented in table 5,2. These data support 
•the view that an ocean color multispectral sensor can do without a 
channel equivalent to channel 1 (0.41 - 0'.46v m) of -the S-192, as no 
information on a yery strong color boundary was contained in those data 
(section 5.3). 

An objective S-192 cloud detection technique was developed 'that- 
uses Gausian statistics to identify cloud-free areas in channel 8 
(0.98 - 1.03 ym) . Cloud-free pixels are 'then analyzed in other channels 
so "tliat contrast stretching based on "the same statistics is automatically 
accomplished (section 5,2). 
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./^pendix A 

This appendix lists the surface data obtained during the 8-9 January 
1974 CGMT) cruises of the R/V VIRGINIA KEY. The geographic positions were 
determined after the control data were carefully replotted and a best fit 
of the vessel"! s location made." Values reported were made using standard 
oceanographic techniques or modified as e^lained in section 2 . 2 of the 
text. 
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TABLE. A.l 
Cl-a 

SURFACE SAMPLES 


3 

a (mg/m } 

TIME (GMT) 

Position 



0.15 

1957 

24° 08.1 


00 

H 

o 

0,09 

2200 

23° 58,8 

N, 

0 
1 — 1 
CO 

0.11 

2400 

23° 47.5' 

N, 

o 

r— 1 
CO 


SAMPLES at depth at 19 57 GMT at 24° 0 8,1 N-, 81° 3 4. 6 W 


• Depth 


Clta (mg/m ) 


10 

m 

0.79 

20 

m 

0,21 

30 

m 

0.25 

40 

m 

0.42 


50 m. 


O'. 30 


SAMPLES at depth at 2200 GMT 
Depth 
10 m 
20 m 
30 m 
40 m 
50 m 


at 23° 58:8 N, 81° 40.4 W 

, 3 

Cl” a Cmg/m, ) 

o-.io • 

0.20 

0.35 

0.31 

0.30 



34. G W 
40.4 W 
47.2 W 



TABLE A. 2 


SALINITY SURFACE SAMPLES 


sal(°/oo) TIME CGMT) Position 


35.996 . 

1314 

24° 

38,9 

N, 

81° 

08.0 

W 

35.924 

1449 

24? 

34.5 

N, 

81° 

12 . 2 

W 

35.929 

1515 

24° 

30.7 

N, 

81° 

16.3 

W 

35.921 

1719 

24° 

23,3 

N, 

81° 

21,8 

W 

36.009 

1806 

24° 

19.1 

N, 

81° 

27.3 

w 

36.023 

1940 

24° 

10.2 

N, 

81° 

33 . 7 

w 

36-. 017 

1958 

24° 

08-.1 

N, 

81° 

34.6 

w 

36,119 

2'200 

23° 

58.8 

N, 

81° 

40,4 

w 

35,874 

2400 

23° 

47.5 

N, 

81° 

47.2 

w 

35.918 

- 0110 

23° 

38.2 

N, 

81° 

52,6 

w 



Table A. 3 ■ 


Volume Scattering .3 (45) 

Surface Samples 


Blue Green GMT 

^ ^ t __ « 


^sr“^ X 10“^) 

(m ^sr"^ X 10”^) 

Time 


Position 



17 . 3 

13.5 

1300 

24*^ 

40.2, 

,N 

81° 

07.3 

W 

13.2 

49.9 

1314 

24° 

38.9 

N 

81° 

08.0 

W 

3,03 

3.02 

1434 

24° 

36.5 

N 

81° 

10 .6 

w 

2.61 

2.08 

1449 

24° 

34.5 

N 

81° 

12.2 

w 

•2,40 

1.60 

1504 

24° 

32.2 

N 

81° 

14.7 

w 

2.82 

2.37 

1515 

24° 

30.7 

N 

81° 

16.3 

w 

2-. 63 

1.95 

1630 

24° 

27.7 

N 

81° 

16.3 

w 

3.07 

2.15 

1645 

24° 

26.1 

N 

81° 

18.1 

w 

.3.. 0.6 

.2 .16 

1700 

24° 

24.6 

N 

81° 

20.2 

w 

2.16 

2.31 

1715 

24° 

23.3 

N 

81° 

21-.8 

w 

3.4.2 

3.09 

1750 

24° 

20.5 

N 

81° 

25.4 

w 

3.10 

2.24 

1802 

24° 

19.1 

N 

81° 

27.3 

w 

4.02 

3.49 

1855 ■ 

24° 

16.3 

N 

81° 

31.0 

w 

2.74 

2.49 

19-10 

24°’ 

14.6 

N 

81° 

32.5 

w 

3.08 

1.75 

1925 

,24° 

12.4 

N 

81° 

33.5 

w 

4.54 ■ 

3.29 

• 1940 

24° 

10.2 

N 

81° 

33.7 

w 

5.82 

8.44 

1957 

24° 

08.1 

N' 

81° 

34.6 

w 

2.76 

2-.- 20 

2040 

24° 

04.5 

N 

81° 

35.6 

w 

2.73 

1.45 

2115 

24° 

03.2 

N 

81°' 

37,7 

w 

3.08 

1.98 

2130 

24° 

01.5 

N 

81° 

39 ,0 

w 

2.99 

• 1.90 

2145 

24° 

00.5 

N 

81° 

39.8 

w 

1.65 

■ 1.77 

2200 

23° 

58.8 

N 

81° 

40.4 

w 


L') 



Table A.* 3 Ccont-’d) 
Volume Sca-ttering 3 <45) 


Blue 

“ Green 

-GMT 


r^sr“^x 10“^) 

(m~^sr“^x 

10“^) Time 


3.35 

2.04 

2255 

23° 

2.63- 

1.85 

2310 

23° 

2.85 

r.83 

2325 

23° 

2'. 58 

1.67 

2400 

23° 

1.92 

' 1.30 

00 50 

23° 

1.93 ■ 

1.28 

0110 

23° 

• 1.96 

1.35 

0130 

23° 

2.13 

1.50 

0150 

23° 


Volume 

Scattering S 

(45) 

Samples- at 

depth at 1314 

GMT at 24° 38 

.9 N, 

depth. 

Blue (m“^sr 

X 10~^) 

Green 


10 m 2.11 


Samples at depth at 1600 GMT at 24° 29.3 N, 
depth Blue x 10~^) _ Green 

10 m 13.. 4 

20m 2.39 

■' 30 m 2.88 

40 m .2.27 

50 m 1.46' 


/.im- 


position 

55.8 N 81° 42.0 W 

53.8 N 81° 43.5 W 
51.3 N 81° 44.6 W 
47.5 N 81° 47.2 W 

40.9 N 81° 51.1 V/ 

38.2 N 81° 52.6 W 
35.7 N 81° 54.1 W 

33.2 N 81° 55.6 W, 


81° 08.0 W 
(m~^sr~^ X 10"“^) , 

I. 83 

81° 15.5 W 

10 "^) 

II . 1 
1.77 
2.80 
1.86 
1,42 



Tab le A. 3 ( cont ’ d ) 
Volume Scattering ^( 4 . 5 ) 


Samples at 

depth 

at 1806 GMT 

at 24° 

19 .1 

N, 

81° 27.3 

W 


depth 

^Blue 

(m“^sr ^ 

X 

10-3) 

Green 

(m-^sr“^ 

X 

10-3) 

10 Jii 


8.07 





8.43 



20 m 


2.57 





2.19 



30 m 


4.37 





3.37 



40 m 


3.34 





2.46 



50 m 


6.87 





5.05 



Samples at 

depth 

at 1957 GMT 

at 24° 

08.1 

N, 

81° 34.6 

W 


depth 

Blue 

(m”lsr*~^ 

X 

10-3) 

Green 

(m“ ^sr- 1 

X 

10-3) 

10 m 


5.68 





3.78 



20 m 


-4.26 





3.66 



30 m 


4.61 





3.73 



40 m 


3.87 





3.26 



50 ‘m 


4.41 





3.55 



Samples at 

depth 

at 2200 GMT 

at 2 3° 

58.8 

N, 

81° 40,4 

W 












depth 

Blue 

(m~^sr ^ 

X 

10-3) 

Green 

(m“ ^sr- ^ 

X 

10-3) 

10 m 


4.12 





2.99 



20 m 


2.60 





1.94 



30 m 


5.70 





4.06 



40 m 


3.98 





3.42 



50- m 


3.28 





2.60 





Temp. (°C) 
24.5 

* 24-..S 

* 25.-3 

* 25.4 

* 25.0 
25.0 

* 24.5 
*24.3 

* 24.4 
*24.5 
*- 24.5 

* 24.7 

* 24.4 
24.4 

* 24.6 
*' 24-. 9 

* 24.9 . 

* 24.8 
25.0 

* 25.0 

* 25.1 

* 25.5 

* 25,5 


Table A. 4 

Bucket Temperatures 


TIME GMT 

Position 






1315 

24° 

38 

.9 

N, 

81° 

08 

.0 

'W 

1434 

24° 

36 

.5 

N, 

81° 

10 

.6 

w 

1449 

24° 

34 

.5 

N, 

81° 

12 

.2 

w 

1504 

24° 

32 

.2 

N, 

81° 

14. 

.7 

w 

1630 

24° 

27 

.7 

M, 

81° 

16 , 

.3 

w 

1645 

24° 

26 

.1 

N, 

81° 

18. 

.1 

w 

1700 

24° 

24 

.6 

N, 

81° 

20 , 

.2 

w 

1750 

24° 

20 

.5 

N, 

81° 

25 . 

.4 

w 

1805 

24° 

19 

.1 

N, 

81° 

27 , 

.3 

w 

1855 

24° 

16 

.3 

N, 

81° 

31, 

.0 

w 

1910 

24° 

14 

.6 

N, 

81° 

32, 

.5 

w 

1925' 

24° 

12 

.4 

N, 

81° 

33. 

.5 

w 

1940 • 

24° 

10 

.2 

N, 

81°- 

33. 

.7 

w 

2010 

24° 

07 

.0 

N, 

81° 

34. 

.6 

w 

2040- 

o 

CM 

04 

.5 

N, 

81° 

35. 

.6. 

. w 

. 2115 

24° 

03 

.2 

N, 

81° 

37. 

.7 

w 

2130 

24° 

01 

.5 

N, 

81° 

39 , 

.0 

w 

2145 

24° 

00 

.5 

N, 

81° 

39. 

.8 

w 

2200 

23° 

58 

.8 

N, 

81° 

40. 

,4 

w 

2255 

23° 

55 

.8 

N, 

81° 

42. 

.0 

w 

2310 

23° 

53 

.8 

N, 

81° 

43, 

,5 

w 

2325 

23°. 

51 

.9 

N, 

81° 

44, 

.6 

w 

2340 

23° 

49 

.9 

N, 

81° 

45, 

,8 

w 

2400 

23° 

47. 

.5 

N, 

81° 

47, 

2 

w 




25.5 




Table A. .4 Ccont.) 
Bucket Temperatures 


Temp. (°C) 

TIME GMT 

Position 





* 26.1 

0050 

23® 

40.9 

N, 

81® 

51.1 

W 

* 26.1 

.0110 

23® 

38.2 

N, 

81® 

52.6 

W 

* 26.2 

0130 

23^ 

' 35.7 

N, 

81® 

54.1 

w 

* 26.3 

0150 

23® 

33.2 

N, 

81®' 

55.6 

w 


* Denotes XBT casts 



APPENDIX B . 


Radiosonde data 

Key West 

1600 GMT, 

8 January 1974 

Pressure (mb) 

Temp. (°C) 

R.H. C%) 


o 



Surface 

1021 

024.8 

79 


-1011 

023.5 

82 

Dry Bulb: 24.°9C 

1000 

023.0 

8'S 

Wet Bulb:' 22.°2C' 

902 

017 .1 

86 

RH: 79% 

-869 

015.9 

63 

Wind Dir: 050°T 

850 

015.1 

64 

Wind Spd: 5 mps 

808 

012.8 

66 

Clouds 1/10 cu • 

799 

012.2 

34 

, 4- 24035 'N 

788 

011.2 

72 ■ 

X 81042’W 

744 

009.9 

66 


740 

008,3 

72 


722 

008.5 

30 


700 

007.8 

31 


671 

005.2 

55 


640 ■ 

003,3 

42 


621 

003-1 

22 


530 

-06.3 

21 


500 

-09,4 

13 


470 

-12.4 

10 


.384 

-24.9 

14 


300 

-39.1 

14 


250 

-48.3 



224 

-50.3 



212 

-48. -9 



200 

-50,3 



150 

-62.6 



100 

-77.1 



070 

-75.7 



066 

-75.2 



061 

-7.1.5 



058- 

-72.7' 



054 

-67.6 



050 

-67.8 



045 

-64.7 



043 

-5'8.1 



• 030 

-51.2 



023 ■ 

-46.3 



020 ■ 

-47.4 



'017 

-48-. 0 



016 

-48.2 



13.5 

-43.9 






APPENDIX C 


This appendix lists the Monte Carlo siinulations of radiances 
and I 2 as described in Section 4.2.2. V7avelengths at which the 
OxN and IxN atinospheric aerosol concentrations were conrputed for. all 
400, 500-, 600, and 78G nm. The cosine of ten zenith angles (p) were 
the independent variables. I, and 1 2 are normalized to unit solar 
flinc on a surface normal to tne solar beam. 



VJave length = M-00 nm 


Aerosol ~ OxN 


V 

.00 

,10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

.95 

1.00 


IlCy) 

.10089+00 

.88850-01 

.78011-01 

.67433-01 

.57839-01 

.52975-01 

.48109-01 

.46607-01 

.44244-01 

.42872-01 

.75043-01 


l2(y) - 

.53255-01 

,72697-01 

.87278-01 

.98883-01 

.10661+00 

.11159+00 

,11406+00 

.116-34+00 

.11615+00 

.11701+00 

.11782+00 


10 



Wavelength = 400 nm 


Aerosol = IxM 


P * 

.00 

,10 

.20 

.30 

.40 

..50 

.60 

.70 

.80 

,90 

.95 

1.00 


IlCy) 

.93753-01 

.90653-01 

.83232-01 

.76449-01 

.67122-01 

.61046-01 

,57118-01 

,54662-01 

.52434-01 

.52278-01 

.71513-01 


I2(P). 

.51310-01 

.63140-01 

.77826-01 

.87640-01 

.96479-01 

.10118+00 

.10457+00 

.10757+00 

.10919+00 

,10916+00 

.10906+00 


II 



Wavelengirh 


500 rim 


Aerosol = OxN 


V 

.00 

.10 

,20 

,30 

.40 

.50 

,60 

.70 


.80 


.90 

.95 

1.00 


.63253-01 

.52971-01 

.37429-01 

.29710-01 

.24246-01 

.22112-01 

.19925-01 

.13015-01 

.18002-01 

.17044-01 

.67308-01 


I2C v) 

.59700-01 

.91821-01 

.10902+00 

.12474+00 

.13138+00 

.13641+00 

,13898+00 

.14171+00 

.14225+00 

,14326+00 

.14203+00 


7 ^ 



Wavelength - 500 


nm 


Aerosol = IxN 


y 

ii(p) 

I 2 C y) 

.00 

. 



.62318-01 

.56038-01 

■ .10 

. 59056-01 

.75227-01 

.20 

.49220-01 

.96058-01 

.30 

.41716-01 

.11022+00 

.40 

.35218-01 

.11781+00 

.50 

. 319287 OI 

.12576+00 

• 

CD 

o 

.27678-01 

,12875+00 

.70 


- 


.26638-01 

.13063+00 

.■80 

.25'5'73-01 

.13228+00 

o 

CD 

. 

.26419-01 

.13315+00 

,95 




.57305-01 

.13155+00 

1C 00 




13 



Wavelength = 600 nm . 


Aerosol = OxN 

a 

U 


12^4) 

.00 

.30110-01 

.42002-01 

.10 

.24508-01 

.80947-01 

.20 

.17292-01 

.10497+00 

' .30 

.13224-01 

.'12254+00 

.40 

.10725-01 

.13037+00 

.50 

■ .98865-02 

.13672+00 

.60 

.90777-02 

.14052+00 

.70 

.88360-02 

.14304+UO 

.80 

.80934-02 

.14358+00 

.90 

.77006-02 

.14664+00 

.95 

1.00 

.61865-01 

.14694+00 


' 7 '/ 



Wavelength = 600 nm 


Aerosol = IxN 


P 

,00 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

,90 

.95 

1,00 


IlCp) 

-.27632-01 

,.3^531-01 

.29628-01 

.25526-01 

.21423-01 

.17841-01 

.16190-01’ 

,15107-01 

.14509-01 

.16199-01 

.52310-01 


I2CP) 

.36991-01 

.67189-01 

.92357-01 

.10823+00 

.11883+00 

.12614+00 

.13059+00 

.13303+00 

.13448+00 

.13524+00 

.13560+00 




780 rm 


Aerosol = OxN 


Wavelength = 


y 

.00 
.10 
.20 
.30 
.40 
.50 
.60 
.70 ■ 
,80 
,90 


2 .95 


.27838-01 

.13026-01 

.76052-02 

.51505-02 

.42167-02 

.37095-02 

.34865-02 

.33889-02 

.30901-02 

.'30811-02 - 

.65909-01 


.69923-01 
.i0854+00 
.13086+00 
.14849+00 
.15502+00 
.16096+00 
.16342+00 
".16269 + 00 
.16464+00 
.16649+00 
.16607+00 


1.00 



Wavelength’ '"= 780 rm 


V 

.00 

.10 

.20 

.30 

.40 

.50 

.60 

r7n 

« « w 
’.'80 


.90 




Aerosol = ! 


l2^v) 

,40506»01 

.67910-01 

.33748*01 

.95346-01 

.24905-01 

.11832+00 

.18224-01 

.13459+00 

.14512-01 

,14180+00 

.11931-01 

-.14856+00 

.10491-01 

.15236+00 

.98166-02 

.15399+00 

.98476-02 

.15395+00 

.11295-01 

.15613+00 

.57149-01 

.15511+00 


1.00 



